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Objectives 


There  has  been  considerable  interest  in  recent  times  in  achieving  consensus  in  swarms  of  agents  using 
simple  strategies.  Work  in  this  area  has  ranged  from  completely  autonomous  agents  to  partially 
autonomous  agents,  with  several  different  frameworks  for  communication  and  information.  In  this 
proposal  we  intended  to  address  the  consensus  problem  where  there  is  a  centralized  decision-maker 
which  conveys  action  plans  to  the  swarm  members  with  the  constraint  that  the  action  plan  is  broadcast  to 
the  members  of  the  swarm,  as  against  individual  control  actions  by  each  agent.  This  implies  that  each 
member  of  the  swarm  receives  exactly  the  same  action  advice  and  implements  the  same.  However,  the 
system  is  not  autonomous  in  the  sense  that  the  information  about  the  agents’  states  is  available  to  the 
centralized  decision-maker  through  its  own  sensing  devices.  We  formulate  the  problem  as  one  in  which 
the  central  decision-maker  has  information  about  the  states  of  the  swarm  members.  This  information  is 
assumed  to  be  obtained  by  the  decision-maker  using  its  own  sensor  suite.  On  the  other  hand,  the  swarm 
members  have  no  global  information.  Their  information  is  local  and  is  restricted  to  a  reduced  state  in 
their  own  local  co-ordinate  system.  They  are  capable  of  receiving  instructions  through  broadcast  and 
implementing  them  in  their  local  reference  frame.  The  problem  we  address  is  to  determine  a  strategy  for 
the  decision-maker  so  that  it  can  broadcast  an  action  plan  to  all  the  swarm  members  at  some  frequency. 
The  objective  is  to  bring  about  a  consensus  among  the  swarm  members  by  having  them  implement  the 
broadcast  instruction.  A  further  objective  was  to  design  a  strategy  that  would  require  minimal 
computational  time. 


Status  of  effort 

In  our  work  we  were  able  to  design  a  strategy  based  on  a  randomization  algorithm  and  a  simplified  linear 
programming  formulation  to  achieve  positional  consensus  among  a  fairly  large  number  of  agents  using  a 
reasonably  small  computational  effort.  Several  strategies  were  examined  and  algorithms  developed  and 
evaluated  and  the  results  have  been  reported.  In  addition  we  have  also  designed  a  consensus  algorithm  in 
which  each  agent  takes  its  decision  autonomously.  We  have  shown  that  the  computational  time 
advantage  is  obtained  even  for  this  case  if  we  resort  to  a  linear  programming  formulation  instead  of  a  the 
standard  non-linear  formulation  that  researchers  have  used  earlier. 


Abstract 

This  work  addresses  the  problem  of  achieving  rendezvous  in  a  multi-agent  system  under  various 
information  paradigms.  We  consider  two  classes  of  algorithms  (i)  Broadcast  based  algorithms  and  (ii) 
Distributed  control  algorithms.  In  the  first  we  consider  both  the  centralized  and  decentralized  cases.  In 
the  centralized  case  each  agent  is  homogeneous  and  all  agents  are  controlled  by  the  same  broadcast 
command  from  a  centralized  controller.  This  method  has  low  communication  cost.  In  the  decentralized 
case  each  agent  computes  its  control  using  the  information  it  obtains  from  its  neighboring  agents  and 
shares  its  control  with  its  neighbours  through  a  broadcast  command.  In  the  second  case  we  consider  each 
agent  to  implement  its  own  control  based  on  information  gathered  from  its  neighbours  through  a  limited 
sensing  capability.  We  show  that  in  the  distributed  control  algorithm,  a  modification  in  the  decision 
domain  of  the  agents  yields  significant  benefits  in  terms  of  computational  time,  when  compared  with 
standard  algorithms  available  in  the  literature.  Moreover,  we  also  show  its  straightforward  application  to 
higher  dimensional  problems  which  is  a  considerable  improvement  over  available  algorithms  in  the 
literature. 

Some  recent  work  in  the  literature,  using  ideal  deterministic  models  for  mobile  robots,  with  its  model 
based  on  an  actual  MEMS  micro-robot,  has  shown  that  it  is  possible  to  achieve  point  convergence  in  this 
framework  for  two  robots  and  to  achieve  limited  consensus.  But  perfect  positional  consensus  cannot  be 
obtained  for  larger  number  of  robots.  In  the  literature,  an  optimization  problem  was  formulated  that 
minimizes  the  maximum  distance  between  swarm  members.  However,  the  formulation  and  the  solution 
suffer  from  two  important  drawbacks.  One  is  that  the  swarm  members  cannot  achieve  point  convergence 
(or  perfect  consensus)  even  after  repeated  application  of  the  strategy.  They  can  only  be  brought  within  a 


certain  distance  of  each  other  and  no  closer.  The  second  drawback  is  that  the  optimization  is  rather  time 
intensive  and  needs  non-standard  optimization  techniques  and  software  such  as  SOCP  (second  order  cone 
programming)  to  obtain  a  solution.  In  our  work  we  propose  a  solution  for  both  the  above  drawbacks.  In 
the  case  of  the  first  drawback  (that  is,  non-achievement  of  consensus)  we  propose  a  strategy  that  has  a 
randomness  function  built  into  the  broadcast  command  that  will  eventually  help  the  swarm  to  achieve 
point  convergence  without  relaxing  the  broadcast  command  paradigm  of  the  original  problem.  In  the 
second  case  (that  is,  large  computation  time),  we  will  propose  a  LP  (linear  programming)  based 
algorithm  that  is  less  computationally  intensive  than  the  SOCP  based  algorithm. 
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Positional  Consensus  in  Multi-Agent  Systems  using  a  Broadcast  Control 

Mechanism 


Kaushik  Das  and  Debasish  Ghose 


Abstract — In  this  paper  a  strategy  for  controlling  a  group 
of  agents  to  achieve  positional  consensus  is  presented.  The 
proposed  technique  is  based  on  the  constraint  that  every  agents 
must  be  given  the  same  control  input  through  a  broadcast 
communication  mechanism.  Although  the  control  command  is 
computed  using  state  information  in  a  global  framework,  the 
control  input  is  implemented  by  the  agents  in  a  local  coordinate 
frame.  We  propose  a  novel  linear  programming  formulation 
that  is  computationally  less  intensive  than  earlier  proposed 
methods.  Moreover,  we  introduce  a  random  perturbation  input 
in  the  control  command  that  helps  us  to  achieve  perfect 
consensus  even  for  a  large  number  of  agents,  which  was  not 
possible  with  the  existing  strategy  in  the  literature.  Moreover, 
we  extend  the  method  to  achieve  positional  consensus  at  a 
pre-specified  location.  The  effectiveness  of  the  approach  is 
illustrated  through  simulation  results. 

I.  Introduction 

The  principle  of  using  multiple  agents  is  motivated  by 
the  idea  that  instead  of  using  a  highly  sophisticated  and 
expensive  robots,  it  may  be  advantageous  in  certain  situations 
to  use  a  group  of  small,  simple,  and  relatively  cheap  robot. 
The  group  of  agents  can  be  used  to  accomplish  various 
tasks  in  different  environment  such  as  tactical  operations, 
exploratory  space  missions,  remote  monitoring  with  mobile 
sensor  networks,  avoidance  of  collision  and  over-crowding 
in  automated  air  traffic  control,  cleanups  of  toxic  spills,  fire 
fighting  and  cooperative  search  with  unmanned  air  vehicles. 

One  of  the  problems  that  is  of  paramount  importance  in 
multi-agent  systems  is  that  of  achieving  consensus,  that  is, 
achieving  identical  values  for  some  specified  subset  of  the 
states  of  the  agents.  For  instance,  the  agents  may  try  to 
converge  to  the  same  direction  of  movement  [1]  after  some 
time  or  they  might  want  to  converge  to  a  point.  Both  are 
problems  in  achieving  consensus.  If  we  have  a  centralized 
system  with  perfect  information  then  achieving  consensus 
is  a  trivial  matter,  since  the  central  controller  can  instruct 
each  agent  suitably  to  reach  a  common  consensus  point. 
However,  if  the  communication  system  has  a  constraint  on 
the  number  of  messages  that  it  can  communicate,  then  one 
may  opt  for  a  broadcast  protocol  where  the  central  controller 
will  communicate  simple  and  identical  instructions  to  all  the 
agents  through  a  broadcast  mechanism.  We  further  impose 
the  additional  constraint  that  each  agent  can  interpret  the 
control  command  only  in  its  local  coordinate  frame  or  local 
state  space.  Only  the  central  controller  has  access  to  the 
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global  states  of  the  system.  Some  of  these  constraint  are 
common  to  other  problems  of  a  similar  nature  (for  instance, 
see  [2],  [3],  [4],  [5]). 

This  problem  was  motivated  by  a  recent  paper  by  Bretl 
[6]  where  a  control  strategy  for  a  group  of  micro-robots  is 
developed  to  perform  a  useful  task  even  when  every  robot 
receives  the  same  control  signal.  The  paper  considers  point 
robots  with  simple  kinematics.  It  was  shown  that  when 
there  are  only  two  agents,  there  exists  a  broadcast  control 
command  (that  is,  both  agents  receive  identical  instructions 
from  the  central  controller)  using  which  both  agents  can  meet 
at  the  same  location  at  the  same  time,  for  almost  all  initial 
conditions.  However,  if  the  number  of  agents  is  more  than 
two,  then  the  best  that  the  agents  can  achieve  is  to  come 
close  to  each  other  within  a  certain  distance  (measured  by 
the  radius  of  the  smallest  disc  that  contains  all  the  agents 
positions),  which  is  a  function  of  the  initial  conditions.  Bretl 
[6]  formulates  this  problem  as  an  optimization  problem  that 
minimizes  the  radius  of  the  disc,  and  proposes  a  solution 
using  the  second  order  cone  programming  (SOCP)  technique 
[8].  However,  using  this  strategy  the  agents  cannot  be  made 
to  converge  to  a  point.  Once  the  solution  of  the  SOCP  is 
implemented,  no  further  improvement  is  possible.  Bretl’s 
paper  was  in  turn  motivated  by  an  interesting  paper  by 
Donald  et  al.  [7]  on  the  development  of  untethered,  steerable 
micro-robots,  where  every  robot  receives  the  same  power  and 
control  signal  through  an  underlying  electrical  grid. 

Our  paper  makes  several  specific  contributions.  The  first 
is  to  propose  a  strategy  that  uses  the  basic  Bretl’s  model 
with  an  additional  randomization  feature  that  allows  large 
number  of  agents  to  achieve  positional  consensus  or  point 
convergence  on  repeated  application  of  the  algorithm  with¬ 
out  compromising  the  broadcast  constraint  on  the  control 
command.  The  second  contribution  is  that  our  method  can 
be  extended  to  the  case  where  the  agents  can  be  made  to 
converge  to  any  pre-specified  point.  The  third  contribution  is 
to  propose  an  optimization  problem  for  this  task  that  is  based 
on  a  linear  programming  formulation.  This  allows  standard 
and  easily  available  software  to  be  used  for  obtaining  the 
solution.  Moreover,  this  formulation  also  retains  the  property 
that  the  number  of  decision  variables,  whose  values  are  to  be 
communicated  to  agents,  remains  unchanged  even  when  the 
number  of  agents  increases.  Finally,  we  also  propose  some 
interesting  properties  related  to  the  positional  formation  of 
agents  when  the  LP  based  strategy  is  applied  iteratively. 

It  is  worth  noting  that  the  randomization  feature  in  the 
algorithm  has  some  similarity  with  the  random  perturbation 
used  in  Viscek’s  model  [1].  Vicsek  et  al.  [1]  propose  a  simple 
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but  compelling  discrete  time  model  of  n  autonomous  agents 
(points  or  particles)  all  moving  in  the  plane  with  constant 
speeds  but  with  different  headings.  Each  agent  updates  its 
heading  using  a  local  rule  based  on  the  average  of  the 
headings  of  its  neighbors  plus  some  random  perturbation. 

The  paper  is  organized  as  follows:  In  Section  II  we 
consider  two  agents  and  show  that  it  is  possible  to  move 
two  agents  to  a  common  location  using  identical  control. 
In  Section  III  we  formulate  a  linear  programming  problem 
for  minimizing  the  proximity  between  agents  by  using  iden¬ 
tical  broadcast  control.  In  Section  IV  we  have  discussed 
some  results  on  the  formation  of  the  agents  after  the  linear 
programming  solution  is  implemented.  In  section  V  we 
introduce  the  notion  of  iterative  solution  of  the  problem  by 
repeated  use  of  the  LP  algorithm  and  show  that  introducing 
a  random  perturbation  in  the  broadcast  mechanism  leads  to 
point  convergence  of  the  agents  by  repeated  application  of 
the  LP  technique.  In  Section  VI  we  present  a  modification 
of  the  algorithm  to  ensure  that  the  swarm  of  agents  converge 
to  a  pre-specified  point.  In  section  VII  we  show  several 
simulation  results  that  illustrate  the  salient  features  of  the 
proposed  algorithm.  Section  VIII  concludes  the  paper  with 
a  discussion  of  possible  future  directions  of  research. 

II.  Formulation  and  Solution  for  Two  Agents 

We  will  first  pose  the  problem  in  a  general  framework 
and  then  address  the  two  agents  case  to  clarify  many  of  the 
assumptions  and  concepts  discussed  in  the  previous  section. 

Assume  that  n  agents  are  located  on  an  obstacle-free 
plane.  We  assume  that  the  central  controller  has  access  to 
the  global  state  of  the  system  which,  in  this  case,  consists 
of  the  position  (x*  £  R2)  and  orientation  ( 9i  £  ( — 7r,  7t] )  of 
the  agents,  i  =  1  The  central  controller  computes 

a  common  local  control  for  the  agents  and  broadcasts  it  to 
the  agents  for  implementation.  The  local  control  is  in  the 
form  of  a  tuple  ( 9,d ),  which  is  interpreted  by  each  agent  in 
its  local  frame  of  reference.  Here,  9  refers  to  the  angle  by 
which  each  agent  changes  its  orientation,  and  d  is  a  scalar 
that  refers  to  the  distance  by  which  each  agent  moves  after 
effecting  the  orientation  change.  Note  that,  the  broadcast 
mechanism  ( 9 ,  d)  is  the  same  for  all  the  agents.  Also,  the 
local  frame  of  reference  for  each  agent  is  centered  at  the 
agent’s  location  and  its  reference  axis  is  oriented  along  its 
current  orientation.  As  an  illustration  see  Fig.  1,  where  agents 
are  shown  located  initially  at  Xjo  with  initial  orientation  9^  in 
the  global  reference  frame.  If  the  control  command  broadcast 
to  all  the  agents  is  ( 9 ,  d),  then  the  agents  implement  it  in 
their  local  coordinate  frame  by  each  of  them  changing  their 
orientation  by  the  same  angle  9  and  advancing  by  the  same 
distance  d  to  reach  the  final  destination  x,f .  Even  in  this 
figure  it  can  be  seen  that  by  doing  this  the  agents  have  come 
closer  to  each  other.  Our  objective  is  to  determine  a  (9,  d) 
such  that  the  agents  can  achieve  the  closest  proximity  with 
each  other. 

Theorem  1:  For  two  agents,  for  all  initial  conditions  of  the 
agents  except  when  9\o  =  6*20,  there  exists  a  control  (9,  d) 
using  which  point  convergence  can  be  achieved. 


Fig.  1.  Basic  configuration 


Note  that  this  result  is  also  available  in  Bretl  [6]  and  is 
stated  here  for  completion.  The  above  theorem  shows  that  it 
is  possible  to  use  a  broadcast  control  command  to  make  two 
agents  meet  at  the  same  location  simultaneously  for  almost 
all  initial  conditions.  However,  the  solution  is  also  unique 
and  hence  the  location  of  the  meeting  point  cannot  be  chosen 
arbitrarily.  One  can  also  interpret  this  result  by  noting  that 
the  final  meeting  point  is  on  the  Voronoi  edge  (equidistant 
line)  between  the  two  initial  positions  of  the  agents.  It  can 
be  shown  that  only  one  unique  point  on  the  Voronoi  edge 
satisfies  the  requirement  that  the  orientation  change  angle 
is  the  same  for  both  the  agents  (see  Fig.  2).  The  point  p 
moves  on  the  equidistant  line  from  —  00  to  +00  and  the 
corresponding  orientation  angle  change  9  is  plotted  for  the 
two  agents.  The  intersection  of  the  two  curves  is  the  unique 
control  command  point.  It  can  be  seen  that  when  the  number 


Fig.  2.  Voronoi  interpretation 


of  agents  is  more  than  two,  each  pair  gives  rise  to  a  different 
unique  meeting  point.  Thus,  there  does  not  exist  a  common 
control  command  to  be  broadcast  so  that  all  the  agents  meet 
at  a  point.  In  the  absence  of  such  a  command,  the  best  that 
can  be  done  is  to  determine  a  (9,  d)  which  brings  the  agents 
in  closest  proximity  with  each  other.  Note  that  in  this  case 
(9,  d)  may  not  be  unique. 

In  the  next  two  sections  we  will  propose  solutions  to 
overcome  both  the  drawbacks  without  compromising  the 
broadcast  based  control  mechanism. 

III.  A  Linear  Programming  Formulation 

Let  the  initial  position  and  initial  orientation  of  the  n 
agents  be  xi0  =  {pa  pi2)  £  R2  and  9t  £  (-7r,7r], 
respectively,  for  all  i  £  {1,  ...,n}.  As  before,  we  define  the 
control  command  to  be  broadcast  as  (0,d).  We  define  our 
performance  measure  as  the  half  length,  denoted  by  r  >  0, 
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of  the  side  of  a  square  oriented  along  the  global  coordinate 
frame,  and  containing  all  the  final  positions  of  the  agents. 
Let  this  square  be  centered  at  2  =  (zi,  Z2)  £  M2. 

Assuming  that  all  the  agents  execute  the  command  (d,d), 
their  final  positions,  given  by  x,/  =  [qn  q *2]  £  R2  will  be, 


That  is, 

qn 

9*2 


Xif 


Xio  +  R(0io)R(0) 


d 

0 


(1) 


Pn 

Pi2 


cos  9io  —  sin  9io 

111 

sin  9io  cos  djo 

.  “2  . 

where,  ui  =  d  cos  9  and  'u2  =  dsind  are  the  control 
variables  that  replace  (9,d).  Note  that  Eqn.  (2)  are  linear 
equations.  Now,  we  formulate  the  linear  programming  prob¬ 
lem  as, 

Minimize  r 
Subject  to 


— r  <  pn  +  ui  cos  0iO  —  u2  sin  9z0  —  z\  <r  (3) 
—r  <  pi2  +  ui  sin  9i0  +  u2  cos  9i0  -  z2  <  r  (4) 
i  =  1, . . . ,  n. 

r>  0  (5) 


of  the  minimal  rectangle  that  contains  the  position  of  all  the 
agents.  Note  that  the  LP  problem  attempts  to  minimize  the 
quantity  r  =  max{5x,  5y}. 

We  first  consider  a  very  special  case  with  three  agents. 
Let  us  assume  that  minimal  formation  by  three  agents  is  a 
rectangle  and  not  necessarily  a  square  where  Sx  =  Sy.  Then, 
there  are  four  way  in  which  this  can  be  occur.  This  is  shown 
in  Figure  3. 


(c)  (d) 


Fig.  3.  Rectangle  formation  by  three  agents  (a)  Two  agents  at  corner,  one 
in  the  interior  (b)  Two  agents  at  corner,  one  on  edge  (c)  All  the  agents  at 
comer  (d)  One  agent  at  a  comer  and  two  agents  on  edges 


The  above  is  a  linear  programming  problem  with  the 
decision  vector  as  (r,  Zi,  z%.  Hi,  112).  Note  that  the  decision 
vector  remains  same  irrespective  of  the  number  of  agents. 
Only  the  number  of  inequality  constraint  increases  with  the 
number  of  agents.  Also,  note  that  Z\,Z2,U\  and  w2  are  free 
variables  and  can  take  both  positive  or  negative  values. 


IV.  Some  Results  on  the  Formation  of  Agents 

After  executing  the  LP  the  distance  between  i  and  j  agent 
along  x-axis  and  y- axis  will  be 


&- 

II 

-  Pji)  +  d(cos(dj  +  9)  -  cos (9j  +  9)) 

=  (Pn 

~  Pji)  +  dC 

(6) 

dy  ij  =  (Pil 

-  Pji)  +  d(sin(9i  +9)  -  sin(#j  +  9)) 

=  (Pil 

l  ~Pji)  +  dS 

(7) 

where,  C  =  (cos(f?j  +  9)  —  cos (Oj  +  9)),  S  =  (sin(&;  +  9)  — 
sin(#j  +9))  and  i.  j  £  {1 . .  .n}.  Thus  the  resulting  distance 
along  x-axis  and  y- axis  are  dictated  by  the  value  of  C  and 
5. 

After  executing  the  LP,  the  new  formation  of  the  agent 
obtained  is  of  interest.  Below,  we  investigate  some  properties 
of  the  formation.  For  this  let  us  define  the  span  of  the 
formation  along  X  and  Y  axis  as  follows:  Let  ( Xi,yi )  be 
the  position  of  the  agents,  where  i  £  I  =  {l,...,n}. 
Then  define  xmax  =  max{xi}iej,  xmin  =  min{xi}i6/, 
Umax  =  max{y,}i6/  and  yrnln  =  min{i/,;},:e/.  Then  the 
span  of  the  formation  along  the  X  and  Y  axis  are  given  by, 

Sx  —  %max  %min 

Sy 

—  Umax  ymin 

The  formation  is  said  to  be  square  if  Sx  =  Sy  and  rectangular 
otherwise.  Essentially,  the  spans  are  the  length  of  the  sides 


Before  we  prove  some  general  results  on  the  formation  of 
the  agents  after  the  LP  is  executed  we  will  state  a  lemma 
that  will  be  useful  to  prove  the  main  results. 

Lemma  1:  If  C  =  cos (9i  +  9)  —  cos (9j  +  9)  then  there 
exists  a  9  such  that  C  <  0  where  9i,9j  £  (— 7r,7r]  and 
9i  7^  9j. 

Proof:  Let,  9  =  (n  —  )  —  A</>).  After  replacing 

0  in  C  =  cos(61;  +  9)  —  cos(9j  +  9)  we  will  get  C  = 
—2  sin(  9‘~9j )  sin(A(ft).  It  is  clear  from  the  expression  that 
we  can  make  C  <  0  by  choosing  A©  properly.  ■ 

Theorem  2:  After  executing  the  LP,  square  is  the  optimal 
formation  for  three  agents. 

Proof:  We  will  prove  this  by  contradiction.  Suppose 
after  executing  the  LP,  the  optimal  formation  is  a  rectangle. 
Let  the  position  and  orientation  of  the  three  agents  be  Xjo  = 
(PiOhPim)  G  R2  and  0,o  £  (— zr, 7r],  respectively,  for  all 
i  £  {1,2,3}.  Let  agents  k  and  l  be  on  the  left  edge  and 
right  edge  of  the  rectangle.  Without  loss  of  generality,  we 
can  assume  Sx 0  >  Sy 0  and  Ad  <  }min{(S,a;o  —  Syo),  \ 
Pioi~Pjoi  | },  where  i,  j  £  (1, 2, 3}.  Let  us  define  a  broadcast 
control  command  (0,  Ad)  such  that  Ad  >  0  is  very  small. 
After  broadcasting  (Ad,  9),  for  some  9,  the  new  positions 
will  be 

Pm  =  Pm  +  Adcos(di0  +  9) 

Pn2  =  Pm2  +  Adsin(#jo  +  9) 

The  new  dimensions  of  the  rectangle  are  given  by  5a;  1  and 
Syi,  where  Sx  1  =  pk  11  -  Pm-  As  Sx0  >  Sy0  and  Ad  < 
\  min{(5x0  -  Sy0),  \  Pioi  -  Pjoi  1}  then  Sy  1  <  Sxl. 

Sx  1  =  Sx0  +  Ad{cos(dfc  +9)  -  cos(9i  +  9)} 

—  Sx  0  +  AdC 
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According  to  Lemma  1,  there  always  exists  9  such  that  C  < 
0.  As  C  <  0  then  5*1  <  Sx o-  This  implies  that  there  exist  a 
broadcast  command  (Ad,  9)  such  that  the  maximum  length 
of  the  sides  of  the  rectangle  can  be  further  reduced.  This 
implies  that  the  optimal  solution  of  the  LP  problem  for  three 
agents  can  not  yield  a  rectangle.  It  has  to  be  a  square. 

This  proof  is  valid  for  Fig.  3(a),  3(b)  and  3(d).  The  proof 
for  Fig.  3(c)  will  be  taken  care  of  by  the  next  Theorem. 

■ 

Theorem  3:  When  the  initial  conditions  are  such  that  only 
one  agent  is  on  one  edge  of  the  longer  span  and  m  agents 
(?n  >1)  are  on  the  other  edge  of  the  larger  span,  then  the 
LP  solution  leads  to  a  square  formation. 


(a)  (b) 

Fig.  4.  An  illustration  for  the  proof  of  Theorem  3. 


Proof:  Consider  Fig.  4  where  Sx  >  Sy.  The  X  distance 
between  agent  j  and  the  other  m  number  of  agent  Sx  = 
Pji—pa,  where  i  £  { 1, . .  .m}.  Let  Ad  <  \(SX  —  Sy).  After 
broadcasting  (Ad,  9),  the  new  dimensions  of  the  rectangle 
are  Sx  and  Sy.  Since  Ad  <  ^(Sx  —  Sy),  we  have  Sy  <  Sx. 
We  can  write 

Sx  =  Pji  -  Pn  +  Ad(cos(9j  —  9)  —  cos (0*  -  9)) 

=  Sx  +  Adf?sin($,;  +  9)  (8) 

where  =  (6j~^et )  and  R  is  a  positive  quantity.  For  a 
particular  03  we  will  get  a  set  of  angles  {<!>;}  for  m  agents. 
Let  $max  =  max {'!»,;}  and  $min  =  ruin  {<!>,}.  Note  that  9, 
is  fixed  here.  Let  9j i  contribute  to  Trnax  and  9j2  contribute 
to  $min.  Then  $max  =  9jl^9i  and  =  9j2+9i .  So 

$max  -  $min  =  9jl~8j2  <  7 t.  We  will  get  Sx  <  Sx,  when 
sin(d>i  +  9)  <0.  The  set  of  angles  {<Jb}  will  be  within  a 
bounded  envelope  of  (0, 7r)  as  4> max  —  $min  £  (0, 7r).  Then 
there  always  exist  an  angle  9  such  that  all  the  envelope  will 
come  in  the  lower  two  quadrants  such  that  all  sin($,;  +  9)  < 
0.  This  implies  that  there  Sx  <  Sx.  Thus  the  LP  solution 
cannot  yield  a  rectangular  formation  since  there  will  be 
another  formation  which  can  be  achieved  by  broadcasting 
and  which  will  have  a  smaller  value  of  max{,S'x;  S'j,}.  ■ 

However,  there  may  not  exist  a  solution  when 
the  number  of  agents  on  the  opposite  edges  are 
more  than  two.  As  an  example  we  can  show 
that  four  agents  with  position  and  orientation  of 
((4.004, 3.9128),  0°),  ((4.015, 2.8264),  5°), ((0.0024, 0.9302) 
,  —  1°)  and  ((0.0049, 1.9007),  0.7°)  will  be  move  to 
((5,4), 5°,  ((5, 3), 10°),  ((1,1),  4°),  ((1,2), 5.7°). 

Although  r  decreases  from  2.0064  to  2  but  Sx  ^  Sy.  The 
reason  behind  this  can  be  explained  as  follows: 


Let  agent  j  and  k  be  on  one  edge  and  a  set  of  agent  { i }  on 
the  opposite  edge  where  i  £  {1, . . .  ,m}.  For  agent  j  there 
will  be  a  set  of  angles  (<I»,;7}  and  for  agent  k  there  will  be 
a  set  of  angles  {4>ifc}.  For  agent  j,  <bmax3  =  max{<FiJ}  and 
T-mm,  =  min{<fbj}  and  for  agent  k,  $maxk  =  max{$ife} 
and  $rnink  =  min{$ifc}.  The  range  of  <l>maXj  -  $minj 
and  $ maxk  -  $mink  will  both  be  (0,7r).  Now,  the  range  of 
max{  }  U  }}  —  min{{$.y  }  U  will  be  greater 

than  (0,7r).  In  which  case  there  does  not  exist  any  common 
9  such  that  the  X  or  Y  axis  distance  will  be  reduced.  They 
will  remain  at  the  same  position  after  executing  the  LP. 

Theorem  4:  If  the  solution  of  the  LP  problem  yields  a 
square  formation  then  the  number  of  agents  on  the  boundary 
of  the  square  is  more  than  two. 

M 


Fig.  5.  An  illustration  for  the  proof  of  Theorem  4. 

Proof:  We  will  prove  this  by  contradiction. Let,  the 
number  of  agent  on  the  square  are  two  and  they  are  located  at 
diagonally  opposite  corners.  The  other  agents  are  the  interior 
of  the  square.  For  the  sake  of  simplicity  we  will  consider 
only  three  agents.  This  is  given  in  Fig.  5, where  MN  is  the 
Voronoi  edge  between  agents  located  at  A  and  B.  According 
to  Theorem  1,  there  always  exists  (9,  d)  such  that  they  can 
meet  at  an  unique  point  on  the  Voronoi  edge  (MN).  Let 
the  unique  meeting  point  be  F.  Let  us  define  a  very  small 
positive  quantity  Ad  such  that  Ad  <  min 772,  ??3, 774}. 
After  broadcasting  (9,  Ad)  the  agents  will  move  from  their 
positions.  The  new  position  of  agents  A  and  B  are  C  and 
D,  respectively.  The  interior  agent  will  remain  in  the  interior 
of  the  square.  CD  is  the  new  diagonal  of  the  square.  We 
can  show  that  A  ABF  ~  A  CDF  and  so  CD\\AB  and 
CD  <  AB.  This  implies  that  further  improvement  of  the 
square  is  possible.  This  is  a  contradiction.  ■ 

V.  Achieving  Perfect  Consensus 

The  solution  of  the  linear  programming  (LP)  problem  will 
yield  control  instructions  which  can  be  broadcast  to  all  the 
agents.  The  agents  will  move  to  a  new  position  or  within  a 
new  square  region  of  smaller  area.  It  can  be  shown  that  no 
further  improvement  of  the  performance  (reduction  in  r)  can 
be  achieved  by  repeated  use  of  the  algorithm.  In  other  words, 
repeated  application  of  the  LP  algorithm  with  the  new  final 
positions  will  not  reduce  the  value  of  r  any  further. 

Suppose  we  represent  the  LP  algorithm  as  an  operation  L 
on  the  initial  conditions  that  yields  the  solution  as, 

L(xi0,9i0\i  =  1  ,...,n)  =  (u\,u*2,r*,Xif,9if)  (9) 
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then, 

L(xif,9if\i  =  1  ,...,n)  =  {0,0,r*,xif,6if)  (10) 

That  is,  there  will  be  no  further  change  in  the  performance 
measure  r.  In  other  words,  ( Xif,9if )  is  a  stationary  point  so 
far  as  the  LP  algorithm  is  concerned. 

We  can  generalize  this  process  by  assuming  that  each  step 
in  the  iteration  is  denoted  by  the  index  k,  with  the  first  step 
in  the  iteration  as  k  =  1.  We  call  this  the  unperturbed  case  as 
the  solution  of  the  LP  is  directly  implemented  by  the  agents 
without  any  perturbation  to  the  solution. 

Theorem  5:  In  the  unperturbed  case,  for  k  >  2,  u\  k  = 
«2  fc  =  0  and  Xi}k+i  =  Xitk',  6(,fc+ i  =  This  means  that 
repeated  use  of  the  LP  solution  on  subsequent  positions  will 
not  reduce  r. 

Proof:  We  will  prove  this  by  contradiction.  Suppose 
for  a  given  initial  condition  (xj0 ,9j0)  we  have  r  =  ro  as 
the  measure  of  proximity  of  the  agents.  Applying  the  LP 
algorithm  we  obtain  u\* ,  v}y ,  rq,  and  the  final  positions 
as  (xh ,  9]j ) .  Now,  considering  the  initial  conditions  as 
(xf0,9f0)  =  {xjf,9jf)  and  applying  the  LP  algorithm  let 
us  assume  that  we  get  itf*  0 7^  0,  r2  <  r±, 
and  the  final  positions  as  (x^,9f^).  Then,  let  us  define 
iti  =  u\*  +  uf*,  z  =  1,2.  It  can  be  shown  that  if  in  the  hrst 
step  Ui  is  used  it  would  yield  a  r  =  r2  <  ?T,  which  implies 
that  r  1  was  not  a  solution  to  the  LP.  This  is  a  contradiction. 

■ 

Now,  consider  a  perturbed  case  where,  the  agents  receive 
a  broadcast  command  containing  the  LP  solution  and  a 
command  to  randomly  perturb  the  final  orientation  angle 
after  the  LP  solution  has  been  implemented.  This  process 
is  shown  in  Figure  6,  where,  the  orientation  angle,  after 
implementing  the  LP  solution  is  perturbed  by  each  agent 
as  follows, 

9i,k+ 1  —  9i^k-\- 1  T  t^i  k+1  (11) 

where,  the  perturbation  angle  Vi^+i  is  given  by  Vi^+i  = 
r/i(3.  rji  is  a  random  number  generated  by  each  agent  inde¬ 
pendently  and  (3  is  a  scaling  angle  which  is  common  to  all 
the  agents.  The  scaling  angle  can  be  set  manually  and  77  can 
be  generated  through  various  distributions.  Here,  we  consider 
both  normal  distribution  and  uniform  distribution. 


Fig.  6.  The  perturbed  case 


VI.  Achieving  Positional  Consensus  at  Desired 
Point 

In  the  previous  section,  we  consider  the  problem  of  posi¬ 
tional  consensus,  but  did  not  have  control  over  at  which  the 


agents  can  meet.  Suppose  we  have  a  pre-specihed  meeting 
point  then  we  can  achieve  this  by  slightly  modifying  the 
previous  formulation.  In  this  modified  form  we  define  the 
meeting  point  as  the  center  of  the  agent  formation  and  is 
denoted  by  (z\,  z2).  Now  the  input  to  the  LP  are  the  initial 
positions,  initial  orientation  and  the  meeting  point.  We  can 
formulate  the  modified  linear  programming  problem  as, 
Minimize  r 
Subject  to 

— t  <  pn  +  Ui  cos  9i0  -  u2  sin  9i0  -  Z\  <  r  (12) 
— r  <  pa  +  u\  sin  9i0  +  u2  cos  9i0  -  z2  <  r  (13) 
i  =  1, . . . ,  n. 

r  >  0  (14) 

The  above  is  a  linear  programming  problem  with  the 
decision  vector  as  (r,ui,u2).  Note  that  the  number  of 
decision  variables  has  reduced  over  the  previous  formula. 
In  the  next  section  we  will  give  simulations  results. 

VII.  Simulation  Results 

In  the  first  set  of  simulations  we  start  with  three 
agents.  Initially  we  consider  ((1, 1),  45°),  ((5, 4),  135°),  and 
((2,6),— 45°)  as  the  initial  position  and  orientation  angle 
of  the  three  agents.  Using  the  perturbation  technique,  with 
normal  distribution  for  77  and  a  scaling  angle  of  j3  =  120°, 
the  agents  converge  to  a  point  after  a  few  iterations  (see 
Fig.  7).  The  variation  in  the  x  and  y  coordinates  of  the 
three  agents  against  the  number  of  iterations  are  also  shown. 
The  convergence  criterion  for  terminating  the  simulation  was 
when  the  value  of  r  became  less  that  2  x  1(L  We  continue 


Fig.  7.  Consensus  with  three  agents  (a)  Trajectory  of  the  agents  (b) 
Reduction  in  r  with  iteration 

with  our  study  with  three  agents.  In  Fig.  8  we  plot  the  average 
number  of  iterations  needed,  and  the  average  length  of  path 
traveled  by  each  agent,  to  achieve  convergence,  as  a  function 
of  the  scaling  angle  (3  for  the  two  cases  when  the  random 
number  77  is  generated  by  a  uniform  distribution  or  by  a 
normal  distribution.  These  results  are  given  by  averaging 
over  200  trials.  We  also  plot  the  standard  deviations.  From 
these  results  we  can  conclude  that  using  larger  scaling  angle 
reduces  r  faster  than  when  the  scaling  angles  are  small.  Also, 
using  normal  distribution  gives  better  results  than  uniform 
distribution. 

Next,  we  consider  larger  number  of  agents  (10  and  15). 
We  use  normal  distribution  and  scaling  angle  of  180°  for 
the  perturbation.  The  convergence  criterion  is  also  relaxed 
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(a)  (b) 


Fig.  8.  Convergence  results  (a)  Average  number  of  iterations  (b)  Average 
length  of  path  (c)  Standard  deviations  for  number  of  iterations:  Uniform 
distribution  (d)  Standard  deviations  for  number  of  iterations:  Normal  dis¬ 
tribution  (e)  Standard  deviations  for  path  length:  Uniform  distribution  (f) 
Standard  deviations  for  path  length:  Normal  distribution 


to  r  <  0.1.  The  results  are  shown  in  Figure  9.  These  results 
show  that  the  computation  time  and  the  number  of  iterations 
rises  with  the  number  of  agents.  This  is  expected  since 
the  complexity  of  the  algorithm  is  same  as  that  of  the  LP 
algorithm. 

To  demonstrate  the  result  that  the  agents  can  be  meet 
any  pre-specified  point,  we  consider  same  set  of  initial 
position  and  orientation  angle(  ((1,1), 45°),  ((5, 4),  135°), 
and  ((2,  6),— 45°)  )  of  the  three  agents.  The  meeting  point 
we  set  as  origin  (0,  0).  The  result  is  illustrated  in  Fig.  10. 

VIII.  Conclusions 

In  this  paper  we  considered  the  problem  of  controlling  a 
group  of  agents  to  converge  at  a  location  using  only  broadcast 
control  input  (identical  control)  for  all  the  agents.  The  results 
shows  that  it  is  possible  for  a  group  of  agents  to  meet  at  a 
location  by  sending  them  a  sequence  of  simple  and  exactly 
identical  instruction.  The  location  point  can  be  pre-specified. 
We  were  able  to  show  that  introducing  a  perturbation  in 
the  broadcast  command  helped  to  achieve  point  convergence 
which  was  not  possible  earlier.  We  also  proposed  a  novel 
linear  programming  based  solution  approach  that  is  compu¬ 
tationally  less  intensive  than  the  SOCP  technique  proposed  in 
the  literature.  There  are  several  opportunities  for  future  work 
in  this  direction.  It  seems  possible  to  extend  this  work  to 
consider  process  noise,  sensor  uncertainty,  and  the  presence 
of  obstacles  in  the  environment.  Moreover,  the  algorithm  can 


Fig.  9.  Convergence  results  (a)  r  vs  number  of  iterations  for  10  agents 
(b)  Trajectory  for  10  agents  (c)  r  vs  number  of  iterations  for  15  agents  (d) 
Trajectory  for  15  agents 


(a)  (b) 

Fig.  10.  Consensus  with  three  agents  at  specified  point  (a)  Trajectory  of 
the  agents  (b)  Reduction  in  r  with  iteration 


most  probably  shown  to  be  robust  to  failures  in  terms  of 
packet  loss  or  failure  of  agents. 
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Abstract:  In  this  paper  we  present  a  modified  broadcast  control  algorithm  for  controlling  a 

group  of  homogeneous  agents  to  achieve  positional  consensus.  We  make  an  assumption  that 
the  agents  are  capable  of  generating  a  flag  which  indicates  whether  agent  will  be  active  or  pas¬ 
sive.  This  proposed  technique  is  based  on  the  constraint  that  all  agents  must  be  given  the  same 
control  input  through  a  broadcast  communication  mechanism.  The  control  command  is  com¬ 
puted  using  state  information  of  active  agents  or  active  and  passive  agents  in  a  global  frame¬ 
work.  The  control  input  is  implemented  by  the  active  agents  in  a  local  coordinate  frame.  We 
propose  some  modified  broadcast  control  mechanisms  based  on  linear  programming  with  or 
without  introducing  random  perturbation  in  the  orientation  of  the  active  agents.  We  show  that 
even  without  introducing  perturbation  in  the  orientation  angle  positional  consensus  is  possible 
if  we  adopt  random  flag  setting.  We  make  a  comparison  between  several  algorithm  that  have 
been  designed  as  combinations  of  the  two  randomization  mechanism.  Moreover,  we  extend 
the  method  to  achieve  positional  consensus  through  this  modified  algorithm  at  a  pre-specified 
point.  The  effectiveness  of  the  approach  is  illustrated  through  simulation  results. 

Key  words:  Multi-agent  systems.  Linear  Programming,  Broadcast  Communications 


1.  INTRODUCTION 

The  principle  of  using  multiple  agents  is  motivated  by  the 
idea  that  instead  of  using  a  highly  sophisticated  and  expen¬ 
sive  agent,  it  may  be  advantageous  in  certain  situations  to  use 
a  group  of  small,  simple,  and  relatively  cheap  agents.  The 
group  of  agents  can  be  used  to  accomplish  various  tasks  in 
different  environment  such  as  tactical  operations,  exploratory 
space  missions,  remote  monitoring  with  mobile  sensor  net¬ 
works,  avoidance  of  collision  and  over-crowding  in  auto¬ 
mated  air  traffic  control,  cleanups  of  toxic  spills,  fire  fighting 
and  cooperative  search  with  unmanned  air  vehicles. 

One  of  the  problems  that  is  of  paramount  importance  in 
multi-agent  systems  is  that  of  achieving  consensus,  that  is, 
achieving  identical  values  for  some  specified  subset  of  the 
states  of  the  agents.  For  instance,  the  agents  may  try  to  con¬ 
verge  to  the  same  direction  of  movement  [1]  or  they  might 
want  to  converge  to  a  point.  Both  are  problems  in  achieving 
consensus.  If  we  have  a  centralized  system  with  perfect  in¬ 
formation,  then  achieving  consensus  is  a  trivial  matter,  since 
the  central  controller  can  instruct  each  agent  suitably  to  reach 
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a  common  consensus  point.  However,  if  the  communication 
system  has  a  constraint  on  the  number  of  messages  that  it 
can  communicate,  then  one  may  opt  for  a  broadcast  proto¬ 
col  where  the  central  controller  will  communicate  simple  and 
identical  instructions  to  all  the  agents  through  a  broadcast 
mechanism.  This  is  the  framework  in  which  we  address  the 
problem  of  consensus  in  this  paper.  Other  papers  that  address 
similar  problem  are  [2,  3,  4,  5], 

This  problem  we  address  in  the  paper  was  motivated  by  a 
recent  paper  by  Bred  [6]  where  a  control  strategy  for  a  group 
of  micro-robots  is  developed  to  perform  a  useful  task  even 
when  every  robot  receives  the  same  control  signal.  The  paper 
considers  point  robots  with  simple  kinematics.  It  was  shown 
that  when  there  are  only  two  agents,  there  exists  a  broad¬ 
cast  control  command  (that  is,  both  agents  receive  identi¬ 
cal  instructions  from  the  central  controller)  using  which  both 
agents  can  meet  at  the  same  location  at  the  same  time,  for  al¬ 
most  all  initial  conditions.  However,  if  the  number  of  agents 
is  more  than  two,  then  the  best  that  the  agents  can  achieve  is 
to  come  close  to  each  other  within  a  certain  distance  (mea¬ 
sured  by  the  radius  of  the  smallest  disc  that  contains  all  the 
agents  positions),  which  is  a  function  of  the  initial  conditions. 
Bretl  [6]  formulates  an  optimization  problem  that  minimizes 
the  radius  of  the  disc,  and  proposes  a  solution  using  the  sec- 
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ond  order  cone  programming  (SOCP)  technique  [7].  How¬ 
ever,  using  this  strategy  the  agents  cannot  be  made  to  con¬ 
verge  to  a  point.  Once  the  solution  of  the  SOCP  is  imple¬ 
mented,  no  further  improvement  is  possible. 

Bretl’s  paper  was  in  turn  motivated  by  an  interesting  pa¬ 
per  by  Donald  et  al.  [8]  on  the  development  of  untethered, 
steerable  micro-robots,  where  every  robot  receives  the  same 
power  and  control  signal  through  an  underlying  electrical 
grid.  These  robots  mainly  have  two  parts.  One  is  an  un¬ 
tethered  scratch  drive  actuator  used  mainly  for  forward  loco¬ 
motion.  The  other  part  is  a  cantilevered  steering  arm  used  for 
turning  through  frictional  contact  with  the  substrate.  These 
micro-robots  are  simple  in  construction  and  can  operate  only 
in  a  local  coordinate  frame.  They  do  not  have  sensory  ca¬ 
pabilities  to  determine  their  position  and  orientation  in  the 
global  coordinate  frame.  Neither  do  they  have  capabilities  to 
localize  themselves  in  relation  to  their  neighboring  objects  or 
other  robots. 

Based  on  the  above  paradigm  Das  and  Ghose  [9]  proposed 
a  strategy  that  uses  the  basic  Bretl’s  model  with  an  addi¬ 
tional  randomization  feature  that  introduced  a  random  per¬ 
turbation  on  the  angular  orientation  of  the  agents  to  achieve 
positional  consensus  or  point  convergence  on  repeated  ap¬ 
plication  of  the  algorithm  without  compromising  the  broad¬ 
cast  constraint  on  the  control  command.  This  method  has 
extended  to  the  case  where  the  agents  can  be  made  to  con¬ 
verge  to  any  pre-specified  point.  The  paper  [9]  also  proposes 
an  optimization  problem  for  this  task  that  is  based  on  a  lin¬ 
ear  programming  formulation  instead  of  SOCP  technique  and 
shows  improved  computational  performance.  However,  the 
proposed  LP  based  method  with  the  randomization  feature 
had  a  serious  drawback  when  the  number  of  agents  is  very 
large.  This  was  observed  in  simulations  carried  out  in  [9].  In 
the  present  paper  we  introduce  another  randomization  feature 
that  not  only  introduces  a  random  perturbation  to  the  agents’ 
angle  but  also  sets  a  flag  that  randomly  assigns  active  or  pas¬ 
sive  role  to  the  agents,  where  an  active  agent  implements  the 
broadcast  control  while  the  passive  agents  ignore  it.  We  show 
that  this  additional  features  helps  in  achieving  consensus  or 
point  convergence  of  even  large  number  of  agents. 


2.  FORMULATION  AND  SOLUTION  FOR  TWO 
AGENTS 


Figure  1:  Basic  configuration 


changes  its  orientation,  and  d  is  a  scalar  that  refers  to  the  dis¬ 
tance  by  which  each  agent  moves  after  effecting  the  orienta¬ 
tion  change.  Note  that  due  to  the  broadcast  mechanism  ( 8 ,  d) 
is  the  same  for  all  the  agents.  Also,  the  local  frame  of  refer¬ 
ence  for  each  agent  is  centered  at  the  agent’s  location  and  its 
reference  axis  is  oriented  along  its  current  orientation.  As  an 
illustration  see  Fig.  1,  where  agents  are  shown  located  ini¬ 
tially  at  Xjo  with  initial  orientation  8,o  in  the  global  reference 
frame. 

If  the  control  command  broadcast  to  all  the  agents  is  (8,  d), 
then  the  agents  implement  it  in  their  local  coordinate  frame 
by  each  of  them  changing  their  orientation  by  the  same  angle 
8  and  advancing  by  the  same  distance  d  to  reach  the  final 
destination  x,/.  Even  in  this  figure  it  can  be  seen  that  by  doing 
this  the  agents  have  come  closer  to  each  other.  Our  objective 
is  to  determine  a  (8,  d)  such  that  the  agents  can  achieve  the 
closest  proximity  with  each  other. 

Theorem  1.  For  two  agents,  for  all  initial  conditions  of  the 
agents  except  when  6\q  =  020,  there  exists  a  control  (8.  d) 
using  which  point  convergence  can  be  achieved. 


Proof.  From  Fig.  1,  assuming  only  two  agents,  we  have. 


xif 


x,o  +  R(8io)R(8 ) 


d 

0 


i  =  1,2 


(1) 


where,  R(a)  is  the  rotation  matrix. 


R(«)  = 


cos  a  -  sin  a 
sin  a  cos  a 


(2) 


We  will  first  pose  the  problem  in  a  general  framework  and 
then  address  the  two  agents  case  to  clarify  many  of  the  as¬ 
sumptions  and  concepts  discussed  in  the  previous  section. 

Assume  that  n  agents  are  located  on  an  obstacle-free  plane. 
We  assume  that  the  central  controller  has  access  to  the  global 
state  of  the  system  which,  in  this  case,  consists  of  the  posi¬ 
tion  (Xj  e  M2)  and  orientation  (8j  e  (-n,  7r])  of  the  agents, 
i  =  1 The  central  controller  computes  a  common  lo¬ 
cal  control  for  the  agents  and  broadcasts  it  to  the  agents  for 
implementation.  The  local  control  is  in  the  form  of  a  tuple 
(8,  d),  which  is  interpreted  by  each  agent  in  its  local  frame 
of  reference.  Here,  8  refers  to  the  angle  by  which  each  agent 


Letting  x\f  -  X2/  we  obtain. 


R(8) 


d 

0 


|7?(02o)  -  R(0io)]  1  (*10  -  x2o) 


(3) 


The  above  equation  can  be  solved  easily  for  d  and  8  if  the 
inverse  on  the  right  hand  side  exists,  which  happens  when 
#10  +  #20-  □ 


The  above  theorem  shows  that  it  is  possible  to  use  a 
broadcast  control  command  to  make  two  agents  meet  at  the 
same  location  simultaneously  for  almost  all  initial  conditions. 
However,  the  solution  is  also  unique  and  hence  the  location 
of  the  meeting  point  cannot  be  chosen  arbitrarily. 
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It  can  be  seen  that  when  the  number  of  agents  is  more  than 
two,  each  pair  gives  rise  to  a  different  unique  meeting  point. 
Thus,  there  does  not  exist  a  common  control  command  to  be 
broadcast  so  that  all  the  agents  meet  at  a  point.  In  the  absence 
of  such  a  command,  the  best  that  can  be  done  is  to  determine 
a  (9,d)  which  brings  the  agents  in  closest  proximity  with  each 
other.  Note  that  in  this  case  (0,  d)  may  not  be  unique. 

In  the  next  two  sections  we  will  propose  a  linear  program¬ 
ming  based  solution  and  a  randomized  mechanism  to  obtain 
consensus  without  compromising  the  broadcast  based  control 
mechanism. 

3.  A  LINEAR  PROGRAMMING  SOLUTION  FOR 
MULTIPLE  AGENTS 

Let  the  initial  position  and  initial  orientation  of  the  n  agents 
be  Xjo  =  (pn  Pa)  e  R2  and  0;  £  ( — zr,  zr],  respectively,  for  all 
i  e  {1 ,...,«}.  As  before,  we  define  the  control  command  to 
be  broadcast  as  (6,  d).  We  define  our  performance  measure 
as  the  half  length,  denoted  by  r  >  0,  of  the  side  of  a  square 
oriented  along  the  global  coordinate  frame,  and  containing  all 
the  final  positions  of  the  agents.  Let  this  square  be  centered 
atz  =  (zi,z2)  e  M2. 

Assuming  that  all  the  agents  execute  the  command  (0,  d), 
their  final  positions,  given  by  x,y  =  [<7,1  qa ]  e  R2  will  be, 

xif  =  x,q  +  R(0jo)R(6)  o  (4) 

That  is, 

<7,i  _  pa  cos  0,-o  -sin  6»,o  u\ 

qa  Pa  sin  6i0  cos  0iO  u2 

where,  ii\  —  d  cos  6  and  u2  —  d  sin  9  are  the  control  variables 
that  replace  (0,  d).  Note  that  Eqn.  (5)  are  linear  equations. 
Now,  we  formulate  the  linear  programming  problem  as. 

Minimize  r 

Subject  to 

—r  <  pa  +  Mi  cos  0/o  -  u2  sin  0,o  -  zi  <  r  (6) 

-r  <  pa  +  mi  sin  0,o  +  u2  cos  0,o  -  z2  <  r  (7) 

i  —  1 

r  >  0  (8) 

The  above  is  a  linear  programming  problem  with  the  deci¬ 
sion  vector  as  (r, zi,Z2,u\,u2).  Note  that  the  decision  vector 
remains  the  same  irrespective  of  the  number  of  agents.  Only 
the  number  of  inequality  constraint  increases  with  the  num¬ 
ber  of  agents.  Also,  note  that  zi,Z2,«i  and  u2  are  free  variables 
and  can  take  both  positive  or  negative  values.  An  illustration 
of  this  process  is  shown  in  Figure  2. 

The  solution  of  the  linear  programming  (LP)  problem  will 
yield  control  instructions  which  can  be  broadcast  to  all  the 
agents  using  which  the  agents  will  move  to  a  new  position  or 
within  a  new  square  region  of  smaller  area,  that  is,  a  square 


Figure  2:  An  illustration  of  how  the  square  reduces  in  size  using  the  LP 
solution 

with  a  reduced  side  length.  It  can  be  shown  that  no  further 
improvement  of  the  performance  (that  is,  reduction  in  r)  can 
be  achieved  by  repeated  use  of  the  algorithm. 

Suppose  we  represent  the  LP  algorithm  as  an  operation  L 
on  the  initial  conditions  that  yields  the  solution  as, 

L(xi0, 0,o  I  *  =  1, . . .  ,n)  =  (u\,  u2,  r*,  xif,  0if)  (9) 

then, 

L{xif,  9jf\i  =  1, . . . ,  n)  =  (0, 0,  r*,  xif,  9if)  (10) 

That  is,  there  will  be  no  further  change  in  the  performance 
measure  r.  In  other  words,  (x,y,  0,/j  is  a  stationary  point  so 
far  as  the  LP  algorithm  is  concerned. 

We  can  generalize  this  process  by  assuming  that  each  step 
in  the  iteration  is  denoted  by  the  index  k,  with  the  first  step  in 
the  iteration  as  k  =  1.  Then  we  can  represent  the  process  as 
in  Figure  3.  We  call  this  the  unperturbed  case  as  the  solution 
of  the  LP  is  directly  implemented  by  the  agents  without  any 
perturbation  to  the  solution.  The  following  result  has  been 
proved  in  [9]  and  is  stated  here  for  completion. 

Theorem  2.  In  the  unperturbed  case,  for  k  >  2,  u\ ,  —  u*2  k  — 
0  and  X/ij+i  =  xty,  0, '_/-+!  =  0,x-  This  means  that  repeated  use 
of  the  LP  solution  on  subsequent  positions  will  not  reduce  r. 


Figure  3:  Generalization  of  the  unperturbed  case 


4.  RANDOMIZATION  MECHANISMS 

Now,  consider  a  perturbed  case  where  the  agents  use  a 
randomization  mechanism  to  modify  their  states  and  imple¬ 
ment  broadcast  command  containing  the  LP  solution.  The 
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broadcast  command  the  new  position  of  the  active  and  pas¬ 
sive  agents  are 


XaJ:+ 1 


Xa,k  +  R(da,k)R(6) 


d 

0 


Xp,k+\  —  %p,k 


where  ( 6 ,  d)  is  the  broadcast  control  command. 


(13) 

(14) 


Figure  4:  Perturbation  in  angle 


5.  ALGORITHMS 


randomization  mechanism  contains  two  components,  in  one 
of  which  a  random  perturbation  is  introduced  in  the  orienta¬ 
tion  angle  of  the  agent  and  in  the  other  a  flag  is  set  by  which 
an  agent  designates  itself  as  active  or  passive.  In  the  angle 
perturbation  case,  the  final  orientation  of  each  agent  is  per¬ 
turbed  and  is  used  by  the  central  controller  to  implement  the 
LP  algorithm.  This  process  is  shown  in  Figure  4,  where,  the 
orientation  angle,  after  implementing  the  LP  solution  is  per¬ 
turbed  by  each  agent  as  follows, 

hk+l  =  Qi.k+ 1  +  Vi  Ml  (11) 

where,  the  perturbation  angle  Vj^+i  is  given  by, 

v,,k+\  =  rliP  (12) 

where,  //,  is  a  random  number  generated  by  each  agent  in¬ 
dependently  and  p  is  a  scaling  angle  which  is  common  to 
all  the  agents.  The  scaling  angle  can  be  set  manually  and  // 
can  be  generated  through  various  distributions.  In  Eqn.  (11) 
v;j,+ 1  represents  the  perturbation  angle.  The  perturbation  an¬ 
gle  for  different  agents  are  different.  This  kind  of  random 
mechanism  has  been  studied  in  [9]  and  has  been  found  to  be 
deficient  when  the  number  of  agents  is  large,  as  they  do  not 
come  closer  than  a  certain  value  even  after  a  large  number  of 
iterations. 

The  main  idea  behind  the  flag  setting  approach  is  that  each 
agent  (i)  randomly  generates  a  flag  j],  that  can  take  one  of 
two  values.  One  is  high  (//,),  denoting  an  active  agent  and 
other  is  low  (//),  denoting  a  passive  agent.  The  set  of  agents 
that  have  high  flags  (active)  is  denoted  by  Ra  and  the  set  of 
agents  that  have  low  flags  (passive)  is  Rp. 

Let  the  set  of  active  agents  at  time  t  be  Ra{t)  =  {i  \  f  =  //,) 
and  the  set  of  passive  agents  at  time  t  be  Rp(t )  =  {;  |  f,  =  //) 
where  i  e  {1 An  active  agent  is  one  that  will  imple¬ 
ment  the  broadcast  command  and  change  its  state  whereas  a 
passive  agent  will  ignore  the  command  and  remain  stationary. 
The  flag  value  for  each  agent  may  or  may  not  be  observed 
by  the  central  controller  (CC).  The  central  controller  have 
information  of  position  and  orientation  of  all  the  agents  (ac¬ 
tive  and  passive)  at  each  step.  The  control  command  ( m*,m *) 
computed  by  the  central  controller  will  be  broadcast  to  all  the 
agents,  but  only  the  active  agents  will  move. 

Let  the  position-orientation  of  active  and  passive  agents 
be  (xfljt,  9a,k)  and  (xPtk,  dP,k)-  respectively.  After  receiving  the 


In  this  section  we  will  propose  five  different  algorithms  ob¬ 
tained  as  combination  of  the  two  randomization  mechanisms 
described  in  the  previous  section.  They  are  denoted  as  A, 
( A,Fo ),  (A.  Lv),  Fo,  Fn  and  are  shown  in  Fig.  5.  The  al¬ 
gorithm  A  that  introduced  perturbation  only  in  the  angle  has 
been  studied  in  [9]. 


CC  :  Central  Controller 


Figure  5:  Taxonomy  of  the  randomization  mechanism  based  algorithm 

5.1.  Algorithm  (A,  Fo) 

In  this  case,  the  active  agents  in  Ra(t)  whose  flag  is  high 
will  perturb  its  orientation  randomly.  Then,  the  central  con¬ 
troller  will  observe  the  positions,  orientations  and  flags  of  all 
agents,  but  will  compute  the  control  command  (6,  d )  using  the 
position  and  orientation  of  the  active  agents  only.  The  control 
command  will  be  broadcast  to  all  the  agents.  After  receiving 
the  control  command  only  the  active  agents  will  change  their 
state.  The  passive  agents  will  remain  stationary.  This  implies 
that, 

LdXatXk))  =  (u\,u*2)  (15) 

where,  0aM  =  0aJ,  +  va,k. 

5.2.  Algorithm  (A,  F^) 

In  this  case  also  the  flag  will  be  set  randomly,  as  well  as 
the  orientation.  However,  the  central  controller  does  not  have 
access  to  the  flag  information  and  will  compute  the  control 
command  (9,  d)  using  positions  and  orientations  of  both  the 
active  agents  and  passive  agents.  The  control  command  will 
be  broadcast  and  all  the  active  agents  will  move  from  their 
respective  positions.  This  implies  that, 

F((ji-a,k->9a,k)i(ji-p,ki9p,k))  —  I  •  ^2 )  (16) 

where,  9ai  =  0aJc  +  va,k. 
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5.3.  Algorithm  ( Fo ) 

In  this  case  the  flags  of  the  agents  will  be  set  randomly 
but  the  orientation  will  not  be  perturbed.  The  central  con¬ 
troller  will  observe  the  flags,  positions  and  orientations  of 
all  the  agents  and  compute  the  control  command  for  the  ac¬ 
tive  agents  only.  The  central  controller  will  broadcast  control 
command  and  all  active  agents  will  move  from  their  respec¬ 
tive  positions.  This  implies  that, 

L(xa,k,6a,k))  =  (u*,u*2)  (17) 

Note  that  in  this  case  there  is  no  perturbation  in  the  agent 
orientation. 

5.4.  Algorithm  ( F a?) 

In  this  case  also  the  flag  of  the  agents  will  be  set  randomly, 
but  orientation  will  not  be  perturbed.  The  central  controller 
will  observe  positions  and  orientations  of  all  the  agents  and 
compute  the  control  command  based  on  the  position  and  ori¬ 
entation  of  all  agents.  The  central  controller  will  broadcast 
the  control  command  to  all  the  agents  but  only  the  active 
agents  will  implement  these  and  will  move. 

T((xfljfc,  @a,k)i  ^p,kT)  —  (u  | ,  r/ 2 1  (18) 

We  put  all  the  above  algorithms  in  a  tabular  form  in  Table  1. 


Table  1 :  The  different  performance  for  the  algorithm  (Y  =  Yes  ;  N  =  No) 


Algorithm 

Random 

Perturba¬ 
tion  in 

Orientation 

Random 

Flag 

Setting 

LP 

Solved 

Using 

Agents 

Move¬ 

ments 

of 

Agents 

A 

Y 

N 

all 

all 

(A,F0) 

Y 

Y 

active 

active 

( A,Fn ) 

Y 

Y 

all 

active 

Fo 

N 

Y 

active 

active 

Fn 

N 

Y 

all 

active 

6.  ACHIEVING  POSITIONAL  CONSENSUS  AT  A  DE¬ 
SIRED  POINT 

In  the  previous  section,  we  considered  the  problem  of  po¬ 
sitional  consensus,  but  did  not  have  control  over  the  point 
at  which  the  agents  can  meet.  In  [9]  a  scheme  has  been  pro¬ 
posed  to  address  this.  We  can  define  the  meeting  point  (z.\ ,  z2) 
of  the  agent  formation.  In  this  case,  the  input  to  the  LP  are 
the  initial  positions,  initial  orientation  and  the  meeting  point. 
The  linear  programming  formulation  is  the  same  as  [9].  We 
can  formulate  the  modified  linear  programming  problem  as, 
Minimize  r 
Subject  to 

-r  <  pai  +  iii  cos  6a0  -  u2  sin  9a0  -  zi  <  r  (19) 

-r  <  pa2  +  ui  sin  9a0  +  u2  cos  9a0  -Zi<r  (20) 

a  =  [i  \  i  €  RJ 

r>  0  (21) 


The  above  is  a  linear  programming  problem  with  the  deci¬ 
sion  vector  as  {r,u\,u2).  Note  that  the  number  of  decision 
variables  has  reduced  over  the  previous  formulation. 

7.  SIMULATION  RESULTS 

In  the  first  set  of  simulations  consider  three  agents. 
Initially  we  consider  ((1,1), 45°),  ((5,4),  135°),  and 

((2, 6),  -45°)  as  the  initial  position  and  orientation  angle  of 
the  three  agents.  Using  algorithm  A  all  three  agents  converge 
to  a  point  after  a  few  iterations  (see  Fig.  6).  The  variation 
in  the  x  and  y  coordinates  of  the  three  agents  against  the 
number  of  iterations  are  also  shown.  The  convergence 
criterion  for  terminating  the  simulation  is  gives  as  value  of  r 
becoming  less  that  2  x  ICC4. 


X-AXIS  ITERATON 


(a)  (b) 


Figure  6:  Consensus  with  three  agents  using  algorithm  A  (a)  Trajectory 
of  the  agents  (b)  Reduction  in  r  with  iteration  (c)  Convergence  in  the  X- 
coordinate  (d)  Convergence  in  the  F-coordinate 

Next,  we  consider  larger  number  of  agents  (5,  10,  15  and 
20).  The  convergence  criterion  is  also  relaxed  to  r  <  0.1. 
We  can  see  after  some  iteration  the  rate  of  decrement  of  per¬ 
formance  radius  r  is  very  small.  For  agents  10,  15,  20  the 
performance  radius  is  not  improving  after  some  iterations. 
The  results  are  shown  in  Figure  7.  These  results  show  that 
the  computation  time  and  the  number  of  iterations  rise  with 
the  number  of  agents.  This  is  expected  since  the  complexity 
of  the  algorithm  is  the  same  as  that  of  the  LP  algorithm. 

Next  we  consider  Algorithm  (A,  F0).  Note  that  in  this  case 
the  control  command  is  computed  by  considering  only  the 
active  agents.  We  consider  5  agents  in  this  simulation.  The 
simulation  result  is  given  in  Fig.  8.  The  convergence  criterion 
for  terminating  the  simulation  is  gives  as  value  of  r  becoming 
less  than  0.1.  The  simulation  result  shows  that  the  agent  are 
not  converging  to  a  point.  In  fact  they  seem  to  be  sometimes 
going  apart  from  each  other  and  sometimes  coming  closer. 
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The  reason  behind  this  is  the  following.  The  central  con¬ 
troller  is  computing  the  control  command  based  on  the  active 
agents  information  only  and  it  is  independent  of  passive  agent 
information.  This  implies  that  the  active  agents  move  to  a 
point  defined  by  the  LP  solution  without  taking  into  consid¬ 
eration  the  presence  of  the  passive  agents.  This  leads  to  the 
active  agents  moving  away  from  the  passive  agents.  So  after 
repeatedly  using  this  algorithm  the  agent  will  not  converge  to 
a  point. 

Next  we  consider  Algorithm  (A,FN).  The  central  con¬ 
troller  does  not  observe  the  flags  and  considers  all  the  agents 
while  compute  the  control  command.  Here  we  consider  5, 10, 
15  and  20  agents  for  illustration.  The  convergence  criterion 
for  terminating  the  simulation  is  the  same  as  (A,  Fo).  The 
simulation  result  is  given  in  Fig.  9.  We  can  see  that  this  al¬ 
gorithm  is  able  to  bring  a  higher  number  of  agents  at  a  point. 
This  algorithm  are  opposite  to  previous  algorithm.  The  rea¬ 
son  why  this  algorithm  works  better  than  (A,  Fo)  is  that  con¬ 
sideration  of  the  passive  agents  in  the  LP  computation  keeps 
the  agents  swarm  together.  Thus,  not  have  access  to  the  flag 
information  (by  the  central  controller)  actually  helps  in  better 
convergence. 

Next  we  consider  the  Algorithm  Fo-  Here  the  central  con¬ 
troller  observes  the  flag  values  and  there  is  no  perturbation 
in  orientation.  Here  also  we  take  5  agents  to  demonstrate 
the  result.  The  termination  condition  for  simulations  is  the 
same.  The  simulation  result  is  given  in  Fig.  10.  We  can  see 
the  agents  do  not  converge  to  a  point.  The  reason  for  not 
converging  is  same  as  Algorithm  (A,Fo)- 

Next  we  consider  the  Algorithm  F^.  Here  the  central  con¬ 
troller  does  not  observes  the  flags  and  there  is  no  perturbation 
of  angles,  we  consider  5,  10,  15  and  20  agents  to  demonstrate 
the  result.  The  simulation  result  are  given  in  Fig.  11.  We  can 
see  the  agents  converge  to  a  point  for  upto  15  agents.  The 
simulation  result  (Fig.  11)  shows  that,  this  algorithm  does 
not  work  for  agent  20  agents  but  still  it  gives  better  result 
than  algorithm  A. 

We  see  all  the  algorithms  are  not  efficient  for  making  the 
agents  converge  to  a  point.  Table  2  shows  performance  com¬ 
parison  of  all  the  algorithms.  Table  3,  Table  4  and  Table  5 
gives  comparison  of  the  number  of  iterations  and  of  the  com¬ 
putations  time  between  the  algorithms  for  10  and  15  and  20 
agents  respectively.  Note  that  NC  denotes  these  cases  where 
the  algorithm  did  not  make  the  agents  converge  to  a  point. 

To  demonstrate  the  result  that  the  agents  can  meet  at  any 
pre-specified  point,  we  consider  5  agents.  Here  the  meeting 
point  is  (0,0)  .  The  result  is  given  in  Fig.  12.  It  shows  that 
where  the  final  point  is  specified  all  the  algorithms  will  work. 
Table  6,  Table  7  and  Table  8  gives  comparison  of  the  number 
of  iterations  and  of  the  computations  time  between  the  al¬ 
gorithms  for  10  and  15  and  20  agents  respectively  when  the 
meeting  point  is  specified.  The  results  shows  that  even  the 
number  of  agents  are  higher,  the  algorithms  (A,  F0 )  and  F0 
are  taking  few  iteration  to  converge. We  consider  50  agents 
also  for  the  algorithms  (A,  Fo),  Fo  and  the  simulation  result 
is  shown  in  Fig.  13. 


(c)  (d) 


No  of  Iteration  X-axis  position 


(e)  (f) 


(g)  (h) 


Figure  7:  Algorithm  A  (a)  r  vs  number  of  iterations  for  5  agents  (b)  Trajec¬ 
tory  for  5  agents  (c)  r  vs  number  of  iterations  for  10  agents  (d)  Trajectory 
for  10  agents  (e)  r  vs  number  of  iterations  for  15  agents  (f)  Trajectory  for 
15  agents  (g)  r  vs  number  of  iterations  for  20  agents  (h)  Trajectory  for  20 
agents 


(a)  (b) 

Figure  8:  Algorithm  (A,  Fo)  (a)  Change  in  r  with  iteration  (b)  Trajectory  of 
the  agents 
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Table  2:  Comparison  between  the  algorithms  for  different  number  of  agents 


Converge 

Algorithm 

5 

10 

15 

20 

A 

Y 

Y 

N 

N 

(A,F0) 

N 

N 

N 

N 

C A,Fn ) 

Y 

Y 

Y 

Y 

Fo 

N 

N 

N 

N 

Fn 

Y 

Y 

Y 

N 

(a)  (b) 


X-axis  position 


Figure  10:  Algorithm  Fo  (a)  Reduction  in  r  with  iteration  (b)  Trajectory  of 
the  agents 


(a)  (b) 


(c) 


•!3o  -5  0  5 

X-axis  position 

(d) 


X-axis  position 


(e) 


(f) 


No  of  Iteration 

(g)  (h) 


X-axis  position 


(a)  (b) 


(c)  (d) 


(e)  (f) 


Figure  9:  Algorithm  (A,  F n)  (a)  r  vs  number  of  iterations  for  5  agents  (b) 
Trajectory  for  5  agents  (c)  r  vs  number  of  iterations  for  10  agents  (d)  Trajec¬ 
tory  for  10  agents  (e)  r  vs  number  of  iterations  for  15  agents  (f)  Trajectory 
for  15  agents  (g)  r  vs  number  of  iterations  for  20  agents  (h)  Trajectory  for 
20  agents 


Figure  1 1 :  Algorithm  Fn  (a)  r  vs  number  of  iterations  for  5  agents  (b)  Tra¬ 
jectory  for  5  agents  (c)  r  vs  number  of  iterations  for  10  agents  (d)  Trajectory 
for  10  agents  (e)  r  vs  number  of  iterations  for  15  agents  (f)  Trajectory  for 
15  agents  (g)  r  vs  number  of  iterations  for  20  agents  (h)  Trajectory  for  20 
agents 
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Table  3:  Comparison  between  algorithm  for  10  agents,  (NC  :  No  Conver¬ 
gence) 


10  Agents 

Algorithm 

Iterations 

Computation  Time 

A 

NC 

NC 

C A,F0 ) 

NC 

NC 

(A,  Fn) 

1262 

26 

Fo 

NC 

NC 

Fn 

1381 

29 

Table  4:  Comparison  between  algorithm  for  15  agents  (NC  :  No  Conver¬ 
gence) 


15  Agents 

Algorithm 

Iterations 

Computation  Time 

A 

NC 

NC 

(A,F0) 

NC 

NC 

C AFn ) 

1874 

59 

Fo 

NC 

NC 

Fn 

5584 

234 

Table  5:  Comparison  between  algorithm  for  20  agents  (NC  :  No  Conver¬ 
gence) 


20  Agents 

Algorithm 

Iterations 

Computation  Time 

A 

NC 

NC 

(A,  F o) 

NC 

NC 

(A,Fn) 

7949 

1056 

Fo 

NC 

NC 

Fn 

NC 

NC 

Table  6:  Comparison  between  algorithm  for  10  agents  with  specified  point 
(NC  :  No  Convergence) 


10  Agents 

Algorithm 

Iterations 

Computation  Time 

A 

NC 

NC 

(A,F0) 

119 

9 

( a,fn ) 

2193 

135 

Fo 

118.25 

8 

Fn 

4101 

229 

Table  7:  Comparison  between  algorithm  for  15  agents  with  specified  point 
(NC  :  No  Convergence) 


15  Agents 

Algorithm 

Iterations 

Computation  Time 

A 

NC 

NC 

(A,  F o) 

235 

10 

C A,Fn ) 

NC 

NC 

Fo 

278 

14 

Fn 

NC 

NC 

Figure  12:  Final  point  is  specified  :  (a)  Algorithm  A  :  r  vs  iteration  (b) 
Algorithm  A  :  trajectory  of  the  agents  (c)  Algorithm  (A,  Fo )  :  r  vs  iteration 
(d)  Algorithm  (A,  F o)  :  trajectory  of  the  agents  (e)  Algorithm  (A,  F^)  :  r  vs 
iteration  (f)  Algorithm  (A,  F n)  :  trajectory  of  the  agents  (g)  Algorithm  Fo  : 
r  vs  iteration  (h)  Algorithm  Fo  :  trajectory  of  the  agents  (i)  Algorithm  F #  : 
r  vs  iteration  (j)  Algorithm  Fpj  :  trajectory  of  the  agents 
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Figure  13:  Final  point  is  specified  :  50  agents  (a)  Algorithm  (A,Fo)  :  r  vs 
iteration  (b)  Algorithm  (A,  F o)  '■  trajectory  of  the  agents  (c)  Algorithm  Fo  : 
r  vs  iteration  (d)  Algorithm  Fq  :  trajectory  of  the  agents 


Table  8:  Comparison  between  algorithm  for  20  agents  with  specified  point 
(NC  :  No  Convergence) 


20  Agents 

Algorithm 

Iterations 

Computation  Time 

A 

NC 

NC 

(A,F0) 

554 

24 

(A,Fn) 

NC 

NC 

Fo 

494 

23 

Fn 

NC 

NC 

8.  CONCLUSIONS 

In  this  paper  we  present  some  modified  broadcast  control 
algorithms.  The  results  shows  that  some  algorithms  are  ef¬ 
fective  for  higher  number  of  agents.  We  have  also  shown  that 
without  perturbation  in  the  orientation,  positional  consensus 
is  possible  by  setting  randomized  flags. 
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Abstract.  The  aim  of  this  paper  is  to  develop  a  computationally  efficient  decen¬ 
tralized  rendezvous  algorithm  for  a  group  of  autonomous  agents.  The  algorithm 
generalizes  the  notion  of  sensor  domain  and  decision  domain  of  agents  to  enable 
implementation  of  simple  computational  algorithms.  Specifically,  the  algorithm 
proposed  in  this  paper  uses  a  rectilinear  decision  domain  (RDD)  as  against  the 
circular  decision  domain  assumed  in  earlier  work.  Because  of  this,  the  compu¬ 
tational  complexity  of  the  algorithm  reduces  considerably  and,  when  compared 
to  the  standard  Ando’s  algorithm  available  in  the  literature,  the  RDD  algorithm 
shows  very  significant  improvement  in  convergence  time  performance.  Analyti¬ 
cal  results  to  prove  convergence  and  supporting  simulation  results  are  presented 
in  the  paper. 

Keywords:  Multi-agent,  Rendezvous,  Decision  domain,  Consensus. 


1  Introduction 


Research  on  multiple  agents  in  the  context  of  robotics  is  motivated  by  the  fact  that 
instead  of  using  a  highly  sophisticated  and  expensive  automated  agent,  it  may  be  ad¬ 
vantageous  to  use  a  group  of  small,  simple,  and  relatively  cheap  autonomous  agents 
(mobile  robots  or  UAVs). 

The  autonomous  control  of  multi  agents  has  emerged  as  a  challenging  problem.  The 
agents  are  assumed  to  have  limited  sensor  and  communication  range  and  execute  some 
local  rule-based  strategy  depending  on  the  information  collected  by  each  agent  from  the 
environment  and  from  neighboring  agents.  One  of  the  generic  tasks  that  such  a  system 
of  agent  is  often  called  upon  to  perform  is  to  physically  bring  all  the  agents  to  a  common 
point.  This  is  called  a  multi-agent  rendezvous  problem  0.  This  problem  is  important 
because  if  rendezvous  is  feasible,  then  more  general  formations  are  also  achievable  El- 
Previous  work  in  this  area  are  by  Ando  et  al.  0  and  Lin  et  al.  0  where  each  agent 
moves  toward  the  rendezvous  point  by  performing  a  sequence  of  “stop-and-go"  moves. 
The  stop  mode  is  basically  the  sensing  period  and  is  an  interval  of  fixed  length.  In  the 
go  mode  the  agents  will  maneuver  in  an  interval  of  variable  length  and  will  move  from 
its  current  position  to  a  new  position.  In  0,  as  also  in  our  paper,  the  sensing  period 
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is  assumed  to  be  zero.  Both  0  and  Si  use  algorithms  that  require  determination  of 
the  smallest  circle  that  contains  a  given  set  of  agent  positions.  These  algorithms  are 
called  “circumcenter  algorithms".  Although  the  complexity  of  this  algorithm  is  proved 
to  be  subexponential  of  order  0(ne<'2+o^'l^2]nn)  0,  the  number  of  actual  computations 
is  fairly  high.  In  our  paper  we  generalize  the  notion  of  sensing  domain  and  decision 
domain  and  show  that  by  using  a  rectilinear  decision  domain  the  computations  can  be 
simplified  considerably,  thus  bringing  down  the  convergence  time. 

We  show  that  our  algorithm  is  far  superior  in  terms  of  computational  time  than 
Ando’s  algorithm  0  which  is  the  standard  algorithm  in  the  literature,  we  consider 
point  robots  with  simple  kinematics  and  instantaneous  directional  motion  as  in  0,  0. 

2  Preliminaries 

Let  R  =  {ai,ci2,  ■  ■  ■  ,a„}  be  the  set  of  robots  or  agents.  The  positions  of  agent  a,-  is 
given  by  /?,  =  (x, , y , )  £  R2.  The  sensor  domain  of  a  agent  a/  is  denoted  as  .S',  and  its 
decision  domain  is  denoted  by  D,  ,  where  Z),  C  .S’,.  Information  sensed  from  the  decision 
domain  is  used  to  implement  the  algorithm.  Essentially,  we  introduce  the  concept  that 
information  from  the  whole  of  the  sensor  domain  need  not  be  used  for  decision-making. 
In  Fig.  |T}i,  we  give  a  schematic  of  these  concepts.  Note  that,  in  general,  p,  need  not  be 
inside  Sj.  An  agent  determines  its  set  of  neighbouring  agents  based  on  D,.  In  this  paper 
we  assume  that  the  sensing  domain  (Si)  of  all  the  agents  is  circular  with  radius  r.  The 
decision  domain  (D,)  is  a  square  of  side  2 d,  with  d  <  -L,  aligned  with  a  pre-specihed 
global  (X .  Y)  reference  frame.  This  is  shown  in  Fig.  |T|i.  The  set  of  neighbors  of  agent 
ai  is  defined  as  Ni  =  {aj  |  (|x,-  —  xy|)  <  d  and  (|y,-  —  yj\)  <  d}. 

Note  that  an  agent  is  also  its  own  neighbor,  so  a,-  £  Ni.  Also,  if  aj  £  Nj  then  a,-  £  Nj. 
In  Ando  et  al.  0  Si  and  Dj  are  the  same  and  are  circles. 


Fig.  1.  (a)  General  sensor  domain  and  decision  domain  (b)  The  sensor  domain  is  circular  and  the 
decision  domain  is  a  square  aligned  with  the  global  X  —  Y  coordinates 


3  Rectilinear  Decision  Domain  (RDD)  Algorithm 

Convergence  to  a  rendezvous  point  in  0  is  proved  through  two  properties:  (i)  Agents 
who  are  neighbours  remain  as  neighbours  (ii)  Agents  come  closer  with  each  other  in 
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Fig.  2.  Agent  a,  will  move  to  the  centroid  of  the  rectangle.  The  square  with  broken  lines 
is  the  decision  domain  of  agent  a,  and  agents  inside  it  belong  to  Nj.  The  solid  rectangle 
is  the  smallest  rectangle  that  contains  the  agents  in  Nj. 


some  sense  till  they  meet  at  a  point.  The  RDD  (rectilinear  decision  domain)  algorithm 
retains  these  two  ideas.  Again  as  in  Q  the  RDD  algorithm  has  the  assumption  that  the 
initial  graph  is  connected. 

Algorithm  RDD  (Rectilinear  Decision  Domain) 

Step  1:  Each  agent  a,  determines  its  neighbour  set  Nj  using  its  decision  domain  D,. 

Step  2:  Each  agent  constructs  the  smallest  rectangle,  aligned  with  the  global  coordinate 
axes,  that  contains  all  the  members  of  its  neighbour  set. 

Step  3:  Each  agent  computes  the  centroid  of  the  rectangle  and  moves  to  it. 

Fig.  ED  illustrates  these  steps  where  Cj(tf)  is  the  centroid  of  the  rectangle  at  the  time 
instant  4.  These  steps  are  similar  to  Ando’s  algorithm  f3|,  but  for  a  few  significant  dif¬ 
ferences.  In  Step  1,  Ando’s  algorithm  determines  neighbours  using  the  sensor  domain 
Sj.  In  Steps  2  and  3,  Ando’s  algorithm  computes  the  circumcenter  of  the  neighbours  and 
moves  toward  it  subject  to  a  constraint.  Unlike  RDD  which  allows  the  agents  to  move 
directly  to  the  centroid,  Anod’s  algorithm  may  not  allow  the  agents  to  reach  the  cen¬ 
troid.  These  two  important  differences  lead  to  high  computational  complexity,  and  thus 
higher  convergence  time,  in  Ando’s  algorithm.  In  RDD,  an  agent  a,  uses  the  informa¬ 
tion  Pj  =  {(xj.yj)  \aj  £  Nj}  where  ctj  £  Nj ,  and  computes  maxjxj},  min{x;},  max{y;} 
and  mill  { y; }  to  obtain  the  rectangle.  The  computational  complexity  of  this  operation  is 
0(n).  We  will  now  state  two  important  theorems. 

Theorem  1.  An  agent’s  movement  will  be  confined  to  a  square  of  side  d  centered  at  the 
agent’s  current  position  and  aligned  with  the  global  reference  frame. 

Proof.  Consider  the  maximum  deviation  of  the  centroid  of  the  rectangular  area  along 
the  X-axis  from  the  agent’s  current  position.  This  will  be  less  than  d/ 2  because  the 
maximum  deviation  of  a  neighbor’s  position  along  the  X-axis  is  d.  Similar  arguments 
hold  for  the  T-axis.  So  the  agent  movement  will  be  confined  within  a  square  of  side  d, 
centered  at  the  agent’s  current  position.  □ 
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Fig.  3.  Agent  a.j,aj  and  its  neighbor  projection  along  A-axis 


Theorem  2.  If  at  any  time  t k,  agents  at  and  aj  are  neighbors,  then  they  will  be  neigh¬ 
bors  for  all  time  t>tk. 

Proof  Let  the  projections  of  agents  a,  and  aj  on  the  X  axis  at  time  4  be  denoted  by 
Xftf)  and  Xj[tf),  respectively,  where  Xj(tk)  >  A,- (4).  Let  the  left  most  and  right  most 
projections  of  the  neighbors  of  agent  a,-  be  Xj  (4)  and  Xf(tk),  respectively.  The  position 

of  agents  a,-  and  aj  on  the  X  axis  at  time  instant  4+i  is  then  Ap(4+i)  =  ;l  A 
where  p  =  i.j.  The  distance  between  the  agents  aj  and  aj  along  the  X  axis  at  time 
instant  4+1  is  |  Xj(tk+i)  -  A,(4+1)  |=  \  \  (X\{tk)  - Xj(tk))  +  (XI(tk) -X[(tk))  |. 

We  can  show  that  |  Xj(tk)  —X](tk)  \<  d,  since  when  X-(tk)  <  xUtk)  <  Xj(tk),  we 
can  write  0  <  Xj(4)  -  Xj (4)  <  A, (4)  -  x!(tk)  <  d.  Again,  when  Xj(4)  <  X\{tk)  < 
Xi(tk),  we  can  write  Xftk)  — X-(tk )  <  0  <  Xj (4 )  —  X-(tk)  <  d.  Similarly,  we  can  prove 
I  X[(tk)  d.  This  implies  that  Xj  (4+1)  —  X,-(4+i)  <  d.  We  can  prove  the  same 

along  the  Y  axis  too.  □ 

The  theorem  proves  that  connectivity  between  agents,  which  plays  an  important  role  in 
the  performance  of  the  algorithm,  is  preserved. 


4  Analysis  of  the  RDD  Algorithm 

Let  the  global  convex  hull  made  by  the  positions  of  the  agents  at  the  time  instant  4  be 
denoted  by  Co(tk).  We  can  define  the  diameter  of  the  convex  hull  at  the  time  instant 
4  as  dia(Co(tk ))  =  max{||  Pi(tk)  —  Pj{tk)  ||},  i,j  £  {1,2,. . .  ,n}.  When  rendezvous  is 
achieved  the  diameter  of  the  global  convex  hull  is  zero.  We  will  first  show  that  the 
diameter  of  the  global  convex  hull  will  reduce  at  each  step. 

Let  us  consider  Coj(tk )  as  the  convex  hull  made  by  the  neighbor  set  of  the  agent  a, 
at  the  time  instant  4.  Let  the  smallest  rectangle  containing  the  neighbor  set  of  agent  a,-, 
and  aligned  along  the  global  coordinate  axes,  at  time  instant  4  be  Rj(tk).  It  is  obvious 
that  Coftf)  C  Rj(tk).  Let  y‘min  =  min{y;},  where  yj  is  the  y— th  coordinate  of  aj  £  Nj. 
Other  variables  are  similarly  defined.  Then  we  have  the  following  result. 

mid(Rj(tk))  =  (1/2 )  ((x1^  +  x'^) ,  (ylmin  +  y‘max)) 


(1) 
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Fig.  4.  Rectangle  made  by  the  neighbor  of  agent  a, 


Theorem  3.  mid  (4,(4))  £  Co,  (4)  and  mid  (4,( 4))  is  not  an  corner  point  ofCoftf} 

Proof.  Consider  Fig.  Q]  Here  4,(4)  is  the  rectangle  ABCD  and  mid  point  of  it  is  M. 
Each  side  must  have  an  agent.  Let  Nab  be  the  set  of  agents  on  side  AB.  Similarly, 
Nbc,  Ncd  and  Nda  are  the  set  of  agents  on  the  side  BC,  CD  and  DA,  respectively.  It 
is  clear  that  Nab  n  Ncd  =  0  and  Nbc  n  Nda  =  0-  Choose  agents  aAB,  aBC,  aCD  and  aAB 
in  such  a  way  that  aAB  £  Nab,  aBC  £  Nbc >  aCD  £  Ncd  and  aDA  £  Nda  ■  The  line  joining 
aAB  and  aBc  will  separate  B  and  M.  Similarly,  the  line  joining  aBc  and  rtcvt  wi  1 1  sep¬ 
arate  B  and  M,  the  line  joining  acD  and  a  da  will  separate  D  and  M  and  line  joining 
oda  and  aAB  will  separate  A  and  M.  So,  we  can  write  M  £  Co(aAB,aBC,aCD,ODA)-  As 
Co(flAB,OBC,ocD,ODA)  C  Co, (4),  then  we  can  write  M  £  Co, (4).  Let,  the  coordinate  of 
the  mid  (4/(4))  is  {xlmid,ylm id).  If  the  rectangle  is  not  a  point  then  it  is  evident  xlmjd  7^ 
4m>  x‘mid  ±  Xlnax >  Ymid  ^  and  fmid  ±  fmax •  This  observation  and  M  £  Co/(4)implies 

that  M  can  not  be  a  corner  point  of  Co,(4).  □ 

Now  consider  any  agent  a,.  Let  the  maximum  distance  along  the  X-axis  on  the  right  side 
between  agent  a,  and  its  neighbors  be  d'rx  and  on  the  left  side  is  d\x.  Similarly,  along  the 
K-axis  the  maximum  distance  above  a,  is  d'ay  and  below  of  a,  is  .The  position  of  the 
agent  a,-  at  time  instant  (4+1)  will  be 

*i(4+i)  =  *i(4)  -  {d\x  -  dlrx)/ 2;  y,(4+i)  =  y,(4)  -  {d‘by  - d‘ay)/ 2  (2) 

The  movement  of  the  agent  a,  will  depends  upon  {d\x,dlrx,dlay,d'hy\ .  The  agent  a,-  will 
be  stationary  if  d\x  =  d'rx  and  dlay  =  d'hy.  Next  we  will  state  a  theorem  that  agents  at  the 
corner  points  global  convex  hull  Co(tf)  cannot  remain  stationary. 

Theorem  4.  For  any  agent  a,-  which  is  at  the  comer  of  Co(tf)  and  has  at  least  one 
non-located  neighbor,  both  d'/x  =  d'rx  and  dlay  =  d'by  cannot  be  satisfied. 

Proof.  We  will  prove  this  by  contradiction.  Let  d\x  =  d'rx  and  d'ay  =  d'hv  for  an  agent  <4 
which  is  at  a  corner  point.  Let  the  agent  a,-  be  at  corner  point  A  in  Lig.  0  The  neighbors 
Uj  and  a  1,  for  which  d\x  =  d'rx  and  d'ay  =  d'by  is  at  B  and  C,  respectively.  Lrom  Lig.0  we 
can  write  d'rx  =  d\  sin  9\ ,  d\x  =  di  sin  62,  dlby  =  d\  cos  0\  and  d‘ay  =  cos  0?-  If  d'rx  =  d\x 
and  d'ay  =  d‘ay. 


di/d2  =  sin  02  /sin  0/  =  cos  02  /  cos  0/ 


(3) 
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Fig.  5.  Agent  a;  is  at  the  corner  point  of  Co{tk) 


Equating  the  above  two  equations  O2  =  0\ .  Thus,  A  will  be  on  the  line  BC.  So,  A  is  not 
a  corner  point.  This  leads  to  the  contradiction.  □ 

From  Theorem  4  and  Theorem  3,  it  is  clear  that  agents  at  the  corner  of  the  convex  will 
move  either  inside  of  the  convex  hull  or  move  along  the  edges  of  the  convex  hull. 

Theorem  5.  No  agent  can  be  at  a  comer  point  of  the  convex  hull  Co(tk)  at  time  4+1. 

Proof  We  prove  this  by  contradiction.  Suppose  an  agent  a,  reaches  the  corner  point 
ci  of  the  convex  hull  Co(tk)  at  time  4+1.  This  implies  that  ci  is  the  mid  point  of  the 
rectangle  made  by  the  neighbors  of  agent  <4  at  time  4.  Then,  according  to  Theorem  3, 
ci  is  inside  the  convex  hull  Co,- (4).  Again  Co, (4)  C  Co(tk).  This  implies  that  ci  is  the 
convex  combination  of  some  points  of  Co(4).  This  leads  to  the  contradiction.  □ 

Theorem  6.  Co  (4+1  )  c  Co(4) 

Proof  Let  there  be  m  corner  points  of  the  convex  hull  Co (4)  given  by  Pc( 4)  = 
{P 1  ,...,pm}  and  m\  corner  points  of  the  convex  hull  Co(4+i)  given  by  Pc(4+i)  = 
{p\,-  ■  . , Pm\  }■  Now,  Pc[tk)  nPc(4+i)  =  0,  because  no  agent  can  reach  at  the  corner  of 
the  convex  hull  Co  (4)  according  to  lemma  5  at  time  instant  4+i.  The  corner  point  agents 
of  Co(4)  must  move  (according  to  Theorem  4).  Since  /?,■  £  Co(tk)  and  /?,-  ^  Co(4+i), 
for  all  pi  £  Pc{tk )•  we  have  Co(4+i)  C  Co(4).  □ 

Next  we  will  show  that  the  diameter  of  the  global  convex  hull  will  reduce  at  each  step. 

Theorem  7.  dia(Co{tk+\)  <  dia(Co(tk)) 

Proof  As  before,  let  Pc(tk)  =  {p\,...,pm}  be  the  set  of  corner  points  of  Co(4)  and 
Pc{tk+  \ )  =  {p\ ,•  •  ■  ,Pmi }  be  the  set  of  corner  points  of  Co(4+i).  According  to  the  defi¬ 
nition  of  the  diameter  of  a  convex  hullfTl 

dia[Co(tf))  =max{\\pi- pj\\,\/pi,pj  £  Pc(tk)}  (4) 

dia(Co{tk+x))  =  max{||p,-p;||,Vp/,p;-  £  ^(4+1)}  (5) 

Again  || Pi~Pj\\  <  dia(Co(tk ))  for  all  p,.pj  £  Co(tk)  and pj,pj  £  Pc{tk).  Here  Pc{tk+X)  £ 
Co(4 )  and  ^(4+1 )  r\Pc(tk)  =  0.  From  this,  we  get  dia{Co{tk+\))  <  dia{Co{tk)).  □ 
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The  above  theorem  tells  us  that  the  sequence  of  convex  hulls,  generated  by  the  positions 
of  the  agents  will  make  a  descending  chain  of  convex  sets.  Under  similar  arguments  as 
Ando  et  al.Q  it  can  be  shown  that  the  convex  hull  will  become  a  point  as  all  the  agents 
converge  to  a  point.  We  omit  the  proof  due  to  space  limitations. 


5  Simulation  Results  and  Implementation  Issues 

In  Fig. [6j  four  snapshots  for  10  agents  have  been  shown.  One  can  see  that  all  the  agents 
eventually  converge  to  a  point.  The  system  converges  when  the  maximum  distance  be¬ 
tween  the  agents  along  the  X  axis  and  Y  axis  is  less  than  the  decision  domain  distance 
(i d ),  since  the  agents  would  converge  to  a  single  point  in  the  very  next  step. 


Fig.  6.  10  agents  placed  randomly  converge  to  a  point 


The  comparative  study  between  Ando’s  algorithm  and  the  RDD  algorithm  is  ex¬ 
ecuted.  In  Table  Q]  the  comparison  of  computational  time  to  converge  is  given.  The 
results  shows  that  RDD  algorithm  is  superior  in  terms  of  computational  time. 

In  Table  ^comparison  of  the  number  of  iterations  to  converge  is  given.  The  iteration 
number  in  case  of  RDD  algorithm  is  a  bit  higher  than  Ando’s  algorithm.  The  reason 
behind  for  this  is  the  decision  domain  in  case  of  RDD  algorithm  is  smaller  than  Ando’s 
algorithm. 

The  algorithm  can  be  implemented  by  any  agent  (robot)  that  can  gather  information 
about  and  from  its  neighbors  (positions  and  orientations)  using  sonar  sensors,  vision 
sensors  or  laser  scanners.  After  generating  the  control  command  the  agents  can  broad¬ 
cast  the  control  command  with  positions  and  orientation  of  the  centroid  information  by 
a  short  range  communication  mechanism.  An  agent  receiving  a  broadcast  control  com¬ 
mand  from  another  agent  which  is  not  in  its  decision  domain  will  ignore  the  command. 
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Table  1.  Comparison  of  Computational  Time 


Time  to  Converge  (sec) 

Number  of 

Ando 

RDD 

Ando/RDD 

Agents 

20 

0.5856s 

0.0329 

17.8065 

50 

3.0677s 

0.1758 

17.4544 

100 

12.1486 

0.5439 

22.3360 

150 

40.7101 

1.7395 

23.4031 

Table  2.  Comparison  of  Number  of  Iterations 


Number  of  Iterations  to  Converge 

Number  of  Ando  RDD  Ando  RDD 


Agents  (Time/Iteration)  (Time/Iteration) 


20 

3.8000 

5.2667 

0.1541 

0.0062 

50 

6.8667 

9.3333 

0.4468 

0.0188 

100 

8.4667 

11.8667 

1.4349 

0.0458 

150 

11.3333 

15.9333 

3.5921 

0.1092 

6  Conclusions 

We  presented  and  analyzed  a  rendezvous  algorithm  considering  a  rectilinear  decision 
domain.  The  computational  complexity  of  RDD  algorithm  is  low  compared  to  the  well 
established  Ando’s  algorithm  in  the  literature.  The  RDD  algorithm  is  simpler  in  terms 
of  few  computation  needed  and  in  relaxing  the  restriction  on  the  movement  of  the  agents 
to  the  centroid  of  the  rectangle.  In  this  sense,  it  is  a  purely  centroidal  algorithm.  The 
simulation  results  support  our  claim. 
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Abstract.  In  this  paper  we  develop  a  Linear  Programming  (LP)  based  decentral¬ 
ized  algorithm  tor  a  group  of  multiple  autonomous  agents  to  achieve  positional 
consensus.  Each  agent  is  capable  of  exchanging  information  about  its  position 
and  orientation  with  other  agents  within  their  sensing  region.  The  method  is  com¬ 
putationally  feasible  and  easy  to  implement.  Analytical  results  are  presented.  The 
etfectiveness  of  the  approach  is  illustrated  with  simulation  results. 


Keywords:  Multi-agent,  Linear  programming.  Consensus,  Decentralize. 


1  Introduction 

In  many  applications  teams  of  robots  or  autonomous  agents  are  being  used  to  achieve  a 
mission  that  was  earlier  entrusted  to  a  single  autonomous  robot.  Such  systems  have  the 
obvious  advantage  of  redundancy,  architectural  simplicity,  and  low  cost.  But  they  are 
mote  difficult  to  control  as  they  have  to  be  controlled  autonomously  through  a  decen¬ 
tralized  control  strategy.  The  main  goal  of  control  strategies  for  groups  of  autonomous 
agents  is  to  achieve  a  coordinated  objective  while  using  local  information.  One  of  the 
problems  that  is  of  great  importance  in  this  field  is  that  of  achieving  consensus,  tha: 
is,  achieving  identical  values  tor  some  specified  subset  of  the  states  of  the  agents.  For 
instance,  the  agents  may  try  to  converge  to  the  same  direction  of  movement  [1]  after 
some  time  or  they  might  want  to  converge  to  a  point.  Both  are  problems  in  achieving 
consensus.  If  we  have  a  centralized  system  with  perfect  information  then  achieving  con¬ 
sensus  is  a  trivial  matter,  since  the  central  controller  can  instruct  each  agent  suitably  t 
reach  a  common  consensus  point.  The  challenging  task  is  to  achieve  consensus  using  a 
distributed  and  decentralized  strategy.  Several  classes  of  problems  in  the  area  has  beer 
addressed  in  the  literature  in  12]  -  [5], 

In  [6]  a  Linear  Programming  ( LP )  based  centralized  formulation  was  proposed  : 
achieve  positional  consensus  of  multiple  agents.  The  central  controller  observes  the  po¬ 
sition  and  orientation  of  the  agents  and,  using  this  information,  it  computes  a  commcr 
contiol  tor  all  the  agents  based  on  a  Linear  Programming  (LP)  technique.  The  centr_ 
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Fig.  1.  (a)  General  sensor  domain  and  decision  domain  (b)  The  sensor  domain  is  circular  and  the 
decision  domain  is  a  square  aligned  with  the  global  X-Y  coordinates 


controller  broadcasts  identical  control  command  to  all  the  agents.  The  agents  imple¬ 
ment  the  common  control  command  in  their  own  local  reference  frame  to  move  toward 
consensus.  It  was  shown  that  perfect  consensus  is  not  possible  by  this  method  when  the 
number  ol  agents  is  more  than  two.  So,  a  post-command  randomization  was  introduced 
for  the  agents.  Using  this  technique  it  was  shown  that  consensus  for  larger  number  of 
agents  could  be  achieved.  However,  the  computation  of  the  command  was  done  in  a 
centralized  framework. 

In  this  paper  we  show  that  it  is  possible  to  achieve  positional  consensus  in  a  de¬ 
centralized  framework,  using  a  similar  LP  based  formulation.  We  assume  the  agents’ 
decision  domain  in  M2  space  agent  is  a  square.  We  assume  that  agents  are  able  to  com¬ 
municate  with  each  other  within  the  sensing  zone  and  they  are  also  able  to  exchange 
information  about  their  position  and  orientation.  Using  these  local  information,  each 
agent  will  run  the  LP  based  algorithm  and  will  compute  a  control  for  itself.  After  com¬ 
puting  the  control  command,  an  agent  can  directly  implement  it  or  it  can  obtain  infor¬ 
mation  from  its  neighbors  about  their  control  commands  and  then  implement  an  average 
control.  We  show  that  the  decentralized  strategy  leads  to  consensus. 


2  Preliminary  Definitions 


Le^  =  ia u «2, ■..,««}  be  the  set  of  robots  or  agents.  The  positions  of  agent  a-,  is 
Xi  ~  (Ptl’Pn)  €  R-  and  the  orientations  is  0,-  G  (-7T,  7r],  The  sensor  domain  of  an  agent 
a,  is  denoted  as  5,  and  its  decision  domain  is  denoted  by  Z9„  where  D,  C  S„  Information 
sensed  from  the  decision  domain  is  used  to  implement  the  algorithm.  Essentially,  we 
introduce  the  concept  that  information  from  the  whole  of  the  sensor  domain  need  not 
be  used  for  decision-making.  In  fact,  we  can  select  decision  domains  that  make  compu¬ 
tations  simpler.  In  Fig.  la,  we  give  a  schematic  of  these  concepts.  Note  that,  in  general, 
P:  need  not  be  inside  S,.  An  agent  determines  its  set  of  neighbouring  agents  based  on 
D,.  In  this  paper  we  assume  that  the  sensing  domain  (S,j  of  all  the  agents  is  circular 
with  radius  r.  The  decision  domain  (£>,-)  is  a  square  of  side  2d,  with  d  <  -J-,  aligned 

with  a  pre-specified  global  (A'.L)  reference  frame.  This  is  shown  in  Fig.  lb.  The  set  of 
neighbors  of  agent  a,  is  defined  as  N(  =  { aj  |  0;-*y|)  <d;  ( |y(-  -  yj |)  <  d ;  i^j}. 
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This  system  can  be  represented  by  a  proximity  graph  G(r)  =  (V.Eit))  where  V 
den°,es  the  „ode  or  vertex  of  ,he  graph.  Each  code  represents  an  agin, «  a  robo,  £(,) 

bSwe^n  th  6  SCt  °  edg?  betWCen  thC  n0dCS  at  any  time'  Let  a‘  ~  ai  denotc  the  edges 
e  agents  a{  and  a,  and  (a,-  aj)  £  E (t)  if  and  only  if  aj  <=  N,.  The  degree  of  any 

node  represents  the  number  of  neighbors  of  the  agents  that  is  associated  with  that  node3! 

3  Algorithm:  A  Linear  Programming  Formulation 

The  LP  based  centralized  consensus  algorithm  is  described  in  [6].  In  which  a  central 
controller  computes  the  control  command  using  the  position  and  orientation  of  the 
agents  and  the  control  is  broadcast  to  all  the  agents.  The  computed  control  will  be  in 

{nffJ  ^t6r  'TCivin8  this  COntro1  command,  each  agent  will  execute  a  rotation 
antole  AQ  and  translation  d.  Here,  we  present  an  algorithm  which  is  decentralized  in 
natui  c. 

1  .et  the  set  of  positions  and  orientation  of  the  neighbors  of  agent  a,  be  {  {Xi  1  f  0  U 

G  Ni'.  10  [6L  We  define  the  contro1  command  computed  by  agent  ^  as 

if  ,h« h,,nn*  “S  nei«hb0l  s  °"'y-  We  our  performance  measure 

as  the  half  length,  denoted  by  r,  >  0,  of  the  side  of  a  square  oriented  along  the  global 

coord, nate  frame,  and  containing  all  the  final  positions  of  the  agent  and  its  neighbors 

,,  th'S  SC|Uare  be  ce"Iered  at  5  =  fe,„Z2,)  e  R2-  Since  all  neighbors  of  agent,,,  execute 
he  same  command  {&(>,. d,).  their  final  positions,  given  hy  xJf  =  [9fl  qa\  €  R2,  are 

(1) 

where  ni,  =  d, cod  16,  ]  and  u2,  =  <f,sm(4a,)  are  the  control  variables  used.  Note  that 
lem  for  agenuf  at  N<W-  W  f°rmula,e  ,he  ,in“r  programming  prob- 


'qa'\ 

Pil 

I”  cos  9jQ  —  sin  0,o 

uh 

<7/2  J 

Pi2  _ 

[  sin  0/o  cos  0/o 

“2 

Minimize 
Subject  to 


~n  <  Pj\  +  u\. cos  Qj  -  U2, sin  6,  - z\,  <  r, 
~ri  <  Pj2  +  u\.  sin  Qj  +  ulj  cos  Q,  -  Z2,  <  r, 
n  >0;  j  e  Nj 


(2) 

(3) 

(4) 


This  is  a  linear  programming  problem  with  the  decision  vector  as  (/■/,*,  .  z, 

Onlv  the  dhC,S1fn  VeCt0!'rCmainS  the  Same  irrcsPective  of  the  number  of  neighbors! 
Only  the  number  of  inequality  constraint  increases  with  the  number  of  agents.  Also 

note  that  Zii,Z2ru\i  and  u2j  are  free  variables  and  can  take  both  +ve  or  -ve. 


4  Decentralized  Control 


In  this  section  we  will  give  a  formal  description  of  the  algorithm.  It  is  obvious  that  the 
utcome  ot  this  algorithm  depends  upon  the  positions  and  orientations  of  the  neighbors 
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of  an  agent  as  well  as  of  its  own.  The  positions  of  the  agents  a,  at  time  instant  /*+1  can 
be  written  as 


*/(**+! )  =  +  U(AGj(tk),dj(tk ))  (5) 

where,  U(A9j,dj)  =  f{x,tk.  Ojtk). 

While  computing  the  control  command,  agent  a ;  may  or  may  not  include  its  own 
position  and  orientation.  Also,  after  computing  the  control  command  it  may  or  may  not 
share  its  control  with  its  neighbor.  This  gives  four  possible  operation  for  implementing 
this  algorithm.  The  four  cases  are  shown  in  the  Table  1 .  When  agent  will  share  control 
command  with  its  neighbors,  we  will  say  that  it  uses  a  broadcast  mechanisms  (which  is 
restricted  to  its  neighbors)  to  communicate  with  its  neighbor.  When  there  is  no  sharing 
of  control  command  then  we  will  say  that  there  is  no  broadcast  communication.  Note 
that  broadcast  based  algorithm  needs  some  effort  in  implementing  the  algorithm  as  the 
broadcast  is  restricted. 


Table  1.  Four  different  cases 


No 

Broadcast 

Broadcast 

Self  Included 
j  E  Nil) {i} 

Case  A 

Case  B 

Self  not  Included 
j  €  N; 

Case  C 

Case  D 

4.1  No  Broadcast:  Case  A  and  Case  C 

The  agents  will  share  its  positions  and  orientation  but  it  will  not  share  its  control  with 
its  neighbor.  It  may  or  may  not  include  its  position  and  orientation  while  computing  the 
control  command.  Under  this  paradigm  there  are  two  cases 

1.  Case  A:  Compute  control  using  the  information  of  the  neighbor  as  well  as  of  its 
own.  Then,  U(AOj'dj)  =  /(*,(/*). 0,-(f*))/e^u{y}. 

2.  Case  C:  Compute  control  using  the  information  of  the  neighbor  only.  U(Adj.dj)  = 

6j(tk))ieNj 

4.2  Broadcasting  -  Averaging  Rule:  Case  B  and  Case  D 

The  agents  will  send  or  broadcast  its  control  to  its  neighbors.  So  a  neighbor  will  have 
more  than  one  control  input.  The  agents  will  compute  an  average  of  all  its  control  and 
will  move  to  a  new  position.  Here  also  two  case  can  arise 

1 .  Case  B:  Compute  the  control  using  information  from  neighbors  as  well  as  of  its  own 

and  broadcast  to  the  neighbors.  Then,  UUl , =^\(Adi(tk),  di(tk))+'ZjzNi(Aej,dj))  • 
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Fig.  2.  An  illusration  of  three  agents 


2.  Case  D:  Compute  control  using  only  the  neighbor  information  and  broadcast  to 


neighbors.  Then,  (/„.  =  — 

ntj 


(AOi(tk),di(tk))  +  Y.jeN,{A0j  .dj) 


5  Modifications:  Preserving  Connectivity 

Applying  the  algorithm  has  the  problem  that  since  the  connectivity  between  agents  may 
not  be  preserved,  as  will  be  seen  from  the  following  example. 

Consider  only  three  agents  which  are  collinearly  placed  at  intervals  of  distance  d,  i.e., 
the  maximum  sensing  distance  (see  Fig.  2)  along  the  X -axis.  Without  loss  of  generality 
we  put  a  boundary  at  d/2  distance  in  the  F-axis.  This  will  ensure  that  the  connectivity 
along  the  F-direction  will  not  be  violated.  From  Fig.  2  we  can  see  that  the  initial 
graph  \s  connected.  The  positions  of  the  agents  at  the  time  instant  t*+i  arc  drdj  and 
a'k.  The  graph  will  not  be  connected  if  6k  >  6j  because  (d/ 2)sin(0*)  <  (d/2)  sin(0y). 
We  can  see  from  the  above  result  that  the  graph  will  not  be  connected  using  the  above 
algorithms  and  positional  consensus  may  not  always  be  possible.  To  ensure  positional 
consensus,  we  have  to  ensure  the  connectivity  of  the  graph  always. 

In  order  to  achieve  connectivity,  the  movement  of  the  agents  should  be  restricted.  We 
call  this  as  pairwise  motion  constraint.  This  is  shown  in  Fig. 3(a).  In  this  figure  agent 
a,  has  one  neighbor  aj.  The  movement  of  agent  a,  is  restricted  within  the  square.  The 
center  of  the  square  is  Qj  which  is  the  mid  point  of  the  line  joining  agents  «,•  and  aj. 
The  side  length  is  d.  This  constraint  ensures  that  if  two  agents  connect  at  time  tk,  then 
they  remain  connected  for  all  time  t  >  tk. 

In  Fig.3(b)  we  can  see  that  agent  «,•  has  two  neighbors  a,  and  ak,  gj  is  the  goal 
point  of  the  agents,  computed  by  the  algorithm.  Now  the  shaded  region  is  the  common 
overlapping  zone  between  the  two  squares.  For  agent  aj,  it  can  move  towards  the  goal 
by  a  distance  dij  and  for  agent  ak  it  can  move  a  distance  dlk.  Eventually,  the  agent  a, 
will  move  the  distance  d  =  min {djj,dik}  towards  the  goal.  So,  in  general, an  agent  will 
move  towards  the  goal  d  —  min {djj},aj  E  Nj. 

The  common  overlap  zone  is  shown  by  the  region  ABCD  in  Fig.4.  The  agent  position 
at  time  instant  tk  is  The  goal  positions  is  at  gi(tk)  which  is  outside  of  the  com¬ 
mon  overlap  zone.  The  agent  can  move  up  to  the  point  af  which  is  at  the  border  of  the 
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Fig.  4.  Agent  movement  in  common  overlap  Region 


common  overlap  region.  To  compute  the  final  positions  of  the  agent «/,  we  have  divided 
the  whole  overlap  zone  into  four  parts. 

1.  In i  :  0g  €  [O,0i)U  [04,2*0 

2.  In2  :  6g  e  [01-02) 

3.  //?3  :  6g  6  [02. 03) 

4.  In4  :  0?  e  [03-04) 


Now  we  define  four  index  function  1 1 , I2 , 13  and  I4.  The  index  7|  can  take  values  from 
{0,1}-  We  can  write 


f  1  ,9g€lni  . 
Xo.e^im  ’ 


/6{1,...,4} 


In  the  four  different  cases  the  movement  of  the  agents  will  be  different.  The  final  posi¬ 
tions  coordinate  of  the  agents  are. 
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Pnf  =  Pn  +  {P2Ci  -  Pa)  * h  * cot 03  +  (P2<3  ~  Pa)  *h* cot 9g 

+(/?2<1  -  Pa)  *  h  *  cot  Q\  +  (p2ci  -  Pa)  *h*  cot  9S  (6) 

P/2,  =  Pa  +  (PiCj  -  P/i )  *  h  *  tan  Qg  +  (pif3  -  p/i )  *  h  *  tan  03 

+(piCl  -  pi i )  *  h  *  tan  9g  +  (piCi  -  p/i )  *  h  *  tan  0i  (7) 

The  above  equations  are  valid  when  the  goal  point  is  outside  common  overlap  region  or 
on  the  boundary  the  common  overlap  region. 

6  Achieving  Perfect  Consensus 

The  solution  of  the  linear  programming  (LP)  based  algorithm  will  yield  control  instruc¬ 
tions  that  can  be  broadcast  to  all  the  neighbors  using  which  the  agents  will  move  to 
a  new  position.  It  can  be  shown  that  no  further  improvement  of  the  performance  (that 
is,  the  agent  will  be  stationary)  can  be  achieved  by  repeated  use  of  the  algorithm.  The 
graph  will  not  change  if  the  set  of  neighbors  do  not  change.  In  other  words,  repeated 
application  of  the  LP  based  decentralized  algorithm  with  the  new  final  positions,  will 
not  yield  any  solution.  Suppose  we  represent  the  LP  algorithm  as  an  operation  L  on  the 
initial  conditions  that  yields  the  solution  as,  L(*/o,0/o|i  =  1, ...,«)  =  {u\.u\,r*,Xjf  ,0//). 

If  the  neighbor  set  remains,  then,  L{x,f.  9jf\i  =  1, . . . ,«)  —  (0.0  ,r*,Xjf.  0,/).  This  was 
the  case  in  the  centralized  strategy  [6].  However,  in  the  decentralized  case,  it  is  possible 
that  the  neighbors  set  of  an  agent  can  change  and  will  yield  non-zero  central  commands 
for  the  agents  at  each  successive  applications  of  the  strategy.  But,  subsequently  when 
agents  come  sufficiently  close  the  graph  becomes  completely  connected  and  the  control 
command  becomes  zero.  In  order  to  continuously  get  non  zero  control  commands,  we 
introduce  a  perturbation  strategy  given  below. 

In  the  perturbed  case  where  the  agents  randomly  perturb  final  orientation  angle  after 
each  LP  solution  is  now  implemented  as  follows  [6], 

0/.jfc+i  =  0/,/t+i  +  v/,*+ 1  (8) 

In  (8)  V,jt+i  represents  the  perturbation  angle. 

7  Simulations  Results  and  Implementation  Issues 

In  the  first  set  of  simulations  we  have  considered  5  agents  which  are  placed  randomly 
at  the  R2  plane.  We  have  set  their  orientation  randomly.  The  result  is  shown  in  Fig.  5. 
We  can  see  that  using  this  algorithm  all  the  agents  eventually  come  to  consensus. 

The  algorithm  can  be  implemented  by  any  agent  (robot)  that  can  gather  information 
about  and  from  its  neighbors  (positions  and  orientations)  using  sonar  sensors,  vision 
sensors  or  laser  scanners.  After  generating  the  control  command  the  agents  can  broad¬ 
cast  the  control  command  with  positions  and  orientation  of  the  centroid  information  by 
a  short  range  communication  mechanism.  An  agent  receiving  a  broadcast  control  com¬ 
mand  from  another  agent  which  is  not  in  its  decision  domain  will  ignore  the  command. 
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Fig.  5.  Modified  Algorithm:  5  agent  preserves  connectivity  at  each  instant  of  time  and  come  closer 
to  a  consensus 


8  Conclusions 

We  presented  a  decentralize  LP  based  algorithm.  This  work  extends  previous  work  on 
broadcast  based  centralized  control  of  agents  to  a  decentralized  framework.  The  flex¬ 
ibility,  that  the  algorithm  gives  to  the  individual  agents  allows  us  to  design  effective 
consensus  algorithms.  In  future,  we  will  focus  on  computational  time  and  implementa¬ 
tion  issues. 
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Abstract 

This  paper  addresses  the  problem  of  achieving  rendezvous  in 
a  multi-agent  system  under  various  information  paradigms. 
We  consider  two  classes  of  algorithms  (i)  Broadcast  based 
algorithms  and  (ii)  Distributed  control  algorithms.  In  the  first 
we  consider  a  case  where  each  agent  is  homogeneous  and  all 
agents  are  controlled  by  the  same  broadcast  command  from  a 
centralized  controller.  This  method  has  low  communication 
cost.  In  the  second  case  we  consider  each  agent  to  implement 
its  own  control  based  on  information  gathered  from  its 
neighbours  through  a  limited  sensing  capability.  In  this  paper 
we  give  a  brief  overview  of  the  broadcast  based  methods  and 
some  results  on  the  distributed  control  algorithm  where  a 
modification  in  the  decision  domain  of  the  agents  is  shown  to 
yield  significant  benefits  in  terms  of  computational  time, 
when  compared  with  standard  algorithms  available  in  the 
literature.  Moreover,  we  also  show  its  straightforward 
application  to  higher  dimensional  problems  which  is  a 
considerable  improvement  over  available  algorithms  in  the 
literature. 

1  Introduction 

Research  on  multiple  agents  is  motivated  by  the  fact  that 
instead  of  using  a  highly  sophisticated  and  expensive 
automated  agent,  it  may  be  advantageous  to  use  a  group  of 
small,  simple,  and  relatively  cheap  autonomous  agents.  The 
group  of  agents  can  be  used  to  accomplish  various  tasks  in 
different  environments  such  as  tactical  operations,  exploratory 
missions,  remote  monitoring  with  mobile  sensor  networks, 
avoidance  of  collision  and  over-crowding  in  automated  air 
traffic  control,  cleanups  of  toxic  spills,  fire  fighting  and 
cooperative  search  with  unmanned  air  vehicles. 

In  multi-agent  systems,  autonomous  control  of  groups  of 
mobile  agents  which  are  loosely  interconnected  through 
limited  range  communication  links  has  emerged  as  a 
challenging  problem.  In  this  paper  we  discussed  two 
information  paradigms;  one  based  on  a  broadcast  mechanism 
and  the  other  on  a  distributed  control  framework.  One  of  the 
generic  tasks  that  such  a  system  of  agent  is  often  called  upon 
to  perform  is  to  physically  bring 


all  the  agents  to  a  common  point.  This  is  called  a  multi-agent 
rendezvous  problem.  This  problem  is  important  because  if 
rendezvous  is  feasible,  then  more  general  formations  are  also 
achievable. 

The  challenge  in  these  algorithms  is  to  develop  control 
algorithms  such  that  the  agents  can  collectively  perform  some 
task.  All  the  agents  have  limited  sensor  and  communication 
range  and  execute  some  local  rule-based  strategy  depending 
on  the  information  collected  by  each  agent  from  the 
environment,  from  neighboring  agents,  and/or  from  a  central 
controller  (in  the  broadcast  case). 

2  Broadcast  Based  Mechanisms 

In  this  framework,  each  agent  executes  a  command  received 
from  a  central  controller  which  computes  a  common 
command  for  all  the  agents  after  observing  their  positions 
and  orientations  [1].  The  main  advantage  of  this  framework 
is  the  saving  in  communication  cost  due  to  the  fact  that  the 
same  command  is  being  transmitted  in  broadcast  mode,  as 
against  individual  commands  to  each  agent.  A  possible 
schematic  representation  is  given  in  Figure  1. 


(a)  (b) 

Figure  1:  (a)  Broadcast  Mechanism  (b)  Implementation  of  the 
broadcast  mechanisms 

Note  that  each  agent  has  a  local  reference  frame  which  is 
used  by  the  agent  to  implement  the  broadcast  command.  For 
instance,  if  the  command  is  “turn  anti-clockwise  by  30 
degrees  and  proceed  for  3  meters”,  then  the  agent  carries  out 
this  command  in  its  own  local  reference  frame  and  not  with 
respect  to  the  global  reference  frame.  Thus,  with  the 


same  command  issued  to  all  agents,  an  agent  that  is  oriented 
at  say  10  deg  from  the  global  reference  frame  will  attain  an 
orientation  angle  of  40  deg,  while  an  agent  with  initial 
orientation  at  1 10  deg  from  the  global  reference  frame  will 
attain  an  orientation  of  140  deg  (Figure  1:).  It  was  shown  in 
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of  agents.  Also,  note  that  r  ,  z,,  It  and  It .  are  free  variables  and  can 
take  both  positive  or  negative  values.  An  illustration  of  this  process  is  shown  in  Figure 
2. 


that  this  problem  can  be  formulated  as  a  linear 
programming  problem,  but  the  best  result  of  bringing  the  agents 
close  can  be  achieved  by  only  one  step.  In  that  paper 

a  random  perturbation  was  introduced  in  each  agent’s 
orientation  and  it  was  shown  that  several  steps  can  then  be  executed 
with  each  step  leading  to  a  smaller  proximity  between  agents  till 
they  converge.  Several  modifications  (based  on  flag  setting  to 
classify  behavior  of  agents)  were  proposed  in  [3]  and  shown  to 
improve  convergence. 

2.1  Linear  Programming  Formulation 


Figure  2:  An  illustration  of  how  the  square  reduces  in 
size  using  the  LP  solution 

2.2  Randomization  Mechanisms 

2.2.1  Perturbation  in  the  orientation 

The  solution  of  the  linear  programming  (LP)  problem  will 
yield  control  instructions  which  can  be  broadcast  to  all  the 
agents  using  which  the  agents  will  move  to  a  new  position  or 
within  a  new  square  region  of  smaller  area.  It  can  be  shown 
that  no  further  improvement  of  the  performance  (that  is, 
reduction  in  r  )  can  be  achieved  by  repeated  use  of  the 
algorithm. 

We  can  generalize  the  process  by  assuming  that  each  step  in 
the  iteration  is  denoted  by  the  index  k  ,  with  the  first  step  in 
the  iteration  as  A:  1 . 

Now,  consider  a  case  where,  the  agents  receive  a  broadcast 
command  containing  the  LP  solution  and  a  command  to 
randomly  perturb  the  final  orientation  angle  after  the  LP  solution 
has  been  implemented.  In  this,  the  orientation  angle,  after 
implementing  the  LP  solution  is  perturbed  by  each  agent  as, 

ilfMitk  ,1 

where,  the  perturbation  angle  |  is  given  by 

ik,l  i 

where,  is  a  random  number  generated  by  each  agent 
independently  and  is  a  scaling  angle  which  is  common  to 
all  the  agents. 

rp  u  sin  m  c$s  z 

i  1, ,  «•  2.2.2  Flag  Setting  Approach 

r  0  The  decision  vector  is  (,rz ,  z  ,u  ,u  ) ,  which  the  same 

v’  l’  2’  l’  2;  ’ 

irrespective  of  the  number  of  agents  remains.  Only  the  The  main  idea  behind  the  flag  setting  approach  is  that  each 

agent  randomly  generates  a  flag  f  that  can  take  one  of 

number  of  inequality  constraint  increases  with  the  number 


Let  the  initial  position  and  initial  orientation  of  the  n 

^  respe^tfv^ly,  for  afl  i  <  {£...,’ «}  .  As  before,  wk 
define  the  control  command  to  be  broadcast  as  (,.  We  define 
our  d  )  performance  measure  as  the 

half  length,  denoted  by  r  0,  of  the  side  of  a  square  oriented 

along  the  global  coordinate  frame,  and  containing  all  the 
final  positions  of  the  agents.  2 

Let  this  square  be  centered  at  z  □  . 

Assuming  that  all  the  agents  execute  the  command  (, ),  d 

]2 

their  final  positions,  given  by  x  IflqU'M  will  be, 

That  is 

?1/Ii8i01  P  sin 

qi{f  m  i'2C0S  i&i()  2 

where,  u  d  cos  and  u  d  sin  are  the  control 

variables  that  replace  (,  ).^lote  that  the  above  equations 
are  linear  .  Now,  we  formulate  the  linear  programming 
problem  as, 

Minimize  r 

Subject  to  rp  u  cos  u  sin  z 

/1 1  102  /Cl 


two  values.  One  is  high 


,  denoting  an  active  agent  and 
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fh 

other  is  low  f  ,  denoting  a  passive  agent.  The  set  of 
agents  that  have  high  flags  (active)  is  denoted  by  1^  and 
agents  that  have  low  flags  (passive)  is  R  . 


Let  the  set  of  active  agents  at  time  t  be  Rt()  {i  \  f  f  ) and 


te^P^ive^ge^s^  \^gre  ,•  >n  An  active 


agent  is  one  that  will  implement  the  broadcast  command  and 
change  its  state  whereas  a  passive  agent  will  ignore  the 
command  and  remain  stationary.  The  flag  value  for  each  agent 
may  or  may  not  be  observed  by  the  central  controller  (CC).  The 
central  controller  has  information  of  position  and  orientation  of 

all  the  agents  (active  and  passive)  at  each 

$  * 

step.  The  control  command  (, ,)  computed  bs/.the  , 

'  12  central  eontroTrer  wilt  be  broadcast  to 

all  the  agents,  but  only  the  active  agents  will  move. 


Let  the  position  and  orientation  of  active  and  passive  agents 
be  ( x  hk  )  jjpd  (x  ,  ) ,  respectively.  After  receiving 

the  broadcast  command  the  new  position  of  the  active  and 
passive  agents  are 

y  'p  X(  )0  ^ 

xak,l  ak,  „ 

XX  U 

pk,  1  pk 

where,  (,  )  &  the  broadcast  control  command. 

Now,  we  will  propose  five  different  algorithms  obtained  as 
combination  of  the  two  randomization  mechanisms  described 
in  the  previous  section.  They  are  denoted  as  A  ,(,  ),(,  J,  F,  Fv 
i IhgtuL  tft.at  introduced  perturbation  only  in 
the  angle  has  been  studied  in  [2]. 


Table  1  :  The  different  performance  for  the  algorithm 
(Y=Yes;  N=No) 


Algorithm 

Random 
perturbation 
in  orientation 

Random 

flag 

setting 

LP  solved 

using 

agents 

Movements 
of  agents 

A 

Y 

N 

All 

All 

(,  o)AF 

Y 

Y 

Active 

Active 

(,  a <)AF 

Y 

Y 

All 

Active 

oF 

N 

Y 

Active 

Active 

nF  2.3  Acli 

N 

ievinp  Pnsitin 

Y 

nal  Conset 

all 

isus  at  Des 

Active 

ired  Point 

either  enters  in  the  sensing  range  of  f 


In  the  previous  section,  we  consider  the  problem  of  positional 
consensus,  but  did  not  have  control  over  the  point  at  which 
the  agents  can  meet. 

Suppose  we  have  a  pre-specified  meeting  point  then  we  can 
achieve  this  by  slightly  modifying  the  previous  formulation. 
In  this  modified  form  we  define  the  meeting  point  as  the 


orientation  and  the  meeting  point.  We  can  formulate  the 
modified  linear  programming  problem  as, 


Minimize  r 

Subject  to  r  p  u  cos  u  sin  z  r 

n  l  /02  mi 

r  pt  u  sin  z  r 

i  1,  ,  n. 
r  0 


CC  :  Central  Controller 


Figure  3:  Taxonomy  of  the  randomization  mechanism  based 
algorithm 


The  decision  vector  w  Lru  m  )  .  Note  that  the  number  of  . 

12  decision  variables  nas  reduced  over  the  previous  formula. 

Distributed  Control  Algorithms 


3  Distributed  Control  Algorithm 

The  present  work  proposes  a  new  modified  algorithm  for 
distributed  implementation  of  a  rendezvous  algorithm  based 
on  a  modified  decision  domain.  Previous  work  in  this  area 
are  done  by  Ando  et  al.  [4]  and  Lin  et  al.  [5]  in  two 
dimensions,  where  each  agent  moves  toward  the  rendezvous 
point  by  performing  a  sequence  of  “stop-and-go”  moves.  The 
stop  mode  is  basically  the  sensing  period  and  is  an  interval  of 
fixed  length.  In  the  go  mode  the  agents  will  maneuver  in  an 
interval  of  variable  length  and  will  move  from  its  current 
position  to  a  new  position.  The  algorithm  proposed  by  Ando 
et  al.  [4]  has  been  extended  to  arbitrary  dimensions  by  Cortes 
et  al.  [6], 
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In  our  paper,  as  in  Ando’s  work,  the  sensing  period  is 
assumed  to  be  zero.  Ando’s  works  use  algorithms  that 
require  determination  of  the  smallest  circle  that  contains  a 
given  set  of  agent  positions.  These  algorithms  are  called 
“circumcenter  algorithms”.  Basically  it  is  a  special  case  with 
dimension  2  of  the  miniball  problem  [7]  in  higher  dimensions 
The  complexity  of  this  problem  is  proved  to  be 

sub-exponential  of  order  ne  in1)  i/41n>wnin 


The  outcome  of  this  algorithm  depends  upon  the  set  of 
neighbors  and  the  relative  distance  between  the  neighbor 
agents.  Hence,  the  position  of  agent  a  at  instant  t  depends 

upon  the  set.of  neighbors,  NtL  j  and  relative  distance 
ik\  between  the  agents  zf(  )  at  rime  instant  f  .bo  we  can 

ii  i 

write  p(t )  =  p(t )  +  U(t ) 

ik  ik  1  ik  1 


d  dimension  of  the  space  with  n  as  the  number  of  agents.  So, 
the  number  of  actual  computations  is  fairly  high.  In  our  work 
we  generalize  the  notion  of  sensing  domain  and  decision 
domain  and  show  that  by  using  a  rectilinear  decision  domain 
the  computations  can  be  simplified  considerably,  thus 
bringing  down  the  convergence  time.  We  show  that  our 
algorithm  is  far  superior  in  terms  of  computational  time  than 
Ando’s  algorithm  which  is  the  standard  algorithm  in  the 
literature.  The  comparison  results  between  RDD  and 
Ando’s  algorithm  in  2D  was  given  in  [4].  In  this  work  we 
extended  it  to  three  dimension  and  the  comparative  results  are 
also  given. 


3.1  Preliminaries 

Let  R^^{,aa,  ,  a  }  be  the  set  of  robots  or  agents.  The 

positions  of  agent  a, is  given  by  p.  (yc.y.)  U,  The 
sensor  domain  of  an  agent  a  is  denoted  as  S  and  its 
decision  domain  is  denoted  by  D  ,  where  D  S .  Information 
sensed  from  the  decision  domain  is  used  to  implement  the 
algorithm.  Essentially,  we  introduce  the  concept  that 
information  from  the  whole  of  the  sensor  domain  need  not  be 
used  for  decision-making.  In  Figure  7  we  give  a  schematic  of 
these  concepts.  Note  that,  in  general,  p  need  not  be  inside  S  . 

An  agent  determines  its  set  of 
neighbouring  agents  based  on  D  .  In  this  work  we  assume 
that  the  sensing  domain  (  S.)  of  all  the  agents  is  circular 
with  radius  r  .  The  decision  domain  D  is  a  square  of 

side  2d ,  with  dr  12  Ji  ,  aligned 

global  (,7file  fere  nee  franmtlThis  is  shown  in  Figure  7. 

The  set  of  neighbors  o f SP VF I ^ ti  n ed  as 

N,  {a  |  (|  x.  x|)  <7and(|>’.  j.|)  d}. 

Note  that  an  agent  is  also  its  own  neighbor,  so  a.  N . 

Also,  if  a  N ’  then  a  N .  In  Ando  et  al.  [5]  S  and 
D  are  the  same  and  are  circles. 

Let  at  t  instant  of  time  the  position  of  agent  a  be 
The (r$lat where  i  ,  n 

position  of  agent  a  with  respect  to  its  neighbor  Cl  (say)  at 
time  instant  t.  is,  where  j  {,i,  } .  m 

ztQ  P  o  Pjkfilk o 


t)  isf^feqPfltfftjtfbliWSbJrfelfiflgg^iifance  and 

k  1 

the  set  of  neighbors.  Hence  we  can  write 


UtQ  fN(t)kzx(ik  Uj‘k  j” 


(a)  (b) 

Figure  4:  (a)  General  sensor  domain  and  decision 
domain  (b)  The  sensor  domain  is  circular  and  the  decision 
domain  is  a  square  aligned  with  the  global 
X  Y  coordinates. 


3.2  Rectilinear  Decision  Domain  (RDD)  Algorithm 

The  RDD  (rectilinear  decision  domain)  algorithm  retains  two 
basic  ideas  for  a  general  rendezvous  problem  (1)  Agents  who 
are  neighbours  remain  as  neighbours  and  (2)  Agents  come 
closer  with  each  other  in  some  sense  till  they  meet  at  a  point. 
The  RDD  algorithm  will  execute  these  following  steps 

1  Each  agent  determines  its  neighbour  set  using  its 
decision  domain. 

2  Each  agent  constructs  the  smallest  rectangle  in 
two  dimensions  (or  cuboids  in  three  dimensions)  aligned 
with  the  global  coordinate  axes  that  contain  all  the 
members  of  its  neighbour  set. 

3  Each  agent  computes  the  centroid  of  the 
rectangle  (or  the  cuboids)  and  moves  to  it. 


To  construct  a  rectangle  or  cuboids  along  the  global  reference 
frame,  an  agent  requires  the  information  of  the  maximum  and 
minimum  coordinates  along  any  axes,  which 
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is  an  one  dimensional  optimization  problem  over  n  points. 

So  the  computational  complexity  of  this  algorithm  is 
n 

Figure  8  illustrates  these  steps  where  ctOis  the  centroid 

/A-orthe  rectangle  at  the  time  instant 

tk.  These  steps  are  similar  to  Ando's  algorithm  [5],  but  for  a 

few  significant  differences.  In  Step  1,  Ando's  algorithm 

determines  neighbours  using  the  sensor  domain  S  .  In  Steps  2 

and  3,  Ando's  algorithm  computes  the  circumcenter  of  the 

neighbours  and  moves  toward  it  subject  to  a  constraint.  Unlike 

RDD  which  allows  the  agents  to  move  directly  to  the  centroid, 

Ando's  algorithm  may  not  allow  the  agents  to  reach  the 

centroid.  These  two  important  differences  lead  to  high 

computational  complexity,  and  thus  higher  convergence  time, 

in  Ando's  algorithm.  In  RDD,  an  agent 

arises  the  information  P  {(.v,v)  au  N} 
where  a  N.,  and  computes  max{  x  },  min{  x .}  ,  ma x{y  } 
and  min{j  }  to  obtain  the  rectangle.  The  computational 
complexity  of  this  operation  is  ()  On 


• 
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Figure  5:  Agents  will  move  to  the  centroid  of  the 
rectangle.  The  square  with  broken  lines  is  the  decision 
domain  of  agent  a  and  agent  inside  it  belong  to  N . 

The  solid  line  rectangle  is  the  smallest  rectangle  that 
contains  the  agent  in  N . 


Agent  a  will  move  to  the  centroid  of  the  rectangle.  The 
square  with  broken  lines  is  the  decision  domain  of  agent  a 
and  agents  inside  it  belong  to  N .  The  solid  rectangle  is  the 
smallest  rectangle  that  contains  the  agents  in  N. . 

It  has  been  proved  in  [4]  that  any  agent  motion  will  be 
confined  within  a  square  of  side  d  centered  at  the  current 
position  and  aligned  with  the  global  reference  frame.  It  is 
also  shown  in  that  paper  [4]  that  if  any  two  agents  are 
neighbors  at  any  time  instant  A  then  they  will  be  neighbor 

for  all  time  instant  ft  . 


3.3  Analysis  of  Rectilinear  Decision  Domain  (RDD) 
Algorithm 


Let  the  global  convex  hull  made  by  the  positions  of  the 

full  I  at  the  time 

instant  tk as  (())  max  □  p  ()  p  (t )} 
dia  Co  t k{  t  U  where, 

ij,  {1,2,  ,«}.  When  rendezvous  is  achieved  the  diameter 

of  the  global  convex  hull  is  zero.  We  will  first  show  that  the 
diameter  of  the  global  convex  hull  will  reduce  at  each  step. 


t©t  a^t&@iJ®tto.4)tull  made  by  the  neighbor  set  of  the  agent  a.at 

the  time  instant  t  .  Let  the  smallest  rectangle  containing  the 

neighbor  set  of  agent  a  ,  and  aligned  along  the  global 

coordinate  axes,  at  time  instant 
'kikiltmQ- It  is  obvious  that  Co  QRt  tR  (t ). 

Let  y  min  { y  }  ,  where  y  is  the  y  coordinate  of  a  N . 

Other  variables  are  similarly  defined.  Then  we  have  the 

following  result. 

mid  R  ((O)  JO) 

It  has  been  proved  in  [4]  that  mid  ((/^ndlm(id) 

(CfflO^fiofca  corner  point  of  t  (). 

Now  consider  any  agent  a  .  Let  the  maximum  distance  along 
the  X  -  axis  on  the  right  side  between  agent  a  and  its 

neighbors  be  d  and  on  the  left  side  is  d  .  Similarly,  i 

along  the  Y  -axis  the  maximum  distance  above  a  is  d  i 

and  below  of  a  is  ft(  .The  position  of  the  agent  a.at  time 
instant  ( t  )  will  be 

xtQ  x  (t )  ( d  d  )I2 , 

ik  ]  ik  lx  rx  Vt ( )  V/(  )  (d  d  )/2 

ik  1  ik  by  ay 

The  movement  of  the  agent  a  will  depend  m  i 
upon  {c/c/ftTWe}agent  a  will  be  stationary  if 

lx  rx  ay  by  ii  ii 

d  d  and  d  d  .  Next  we  will  use  a  result  that 

Ixrx  ayby 

Co  (.agents  at  the  corner  points  of  global  convex  hull  () 
cannot  remain  stationary. 

It  can  be  proved  that  for  any  agent  a  which  is  at  the  comer  of 
()cfrfdte£  at  least  one  non-located  neighbor,  both 
d  d  and  d  d  cannot  be  satisfied. 

Ixrx  ayby 

Now,  it  is  clear  that  agents  at  the  corner  of  the  convex  will 
move  either  inside  of  the  convex  hull  or  move  along  the 
edges  of  the  convex  hull.  It  has  been  proved  in  [4]  that  Co  t( 
)  Co  t  (  t) .  It  is  shown  in  [4]  that  the  diameter 


ISBN  978-E 


of  the  global  convex  hull  will  reduce  at  each  step,  so  (( t  ^ 
dia  Co  t  (j)The  above  result  tells  us  that  the  sequence  of 
convex  hulls,  generated  by  the  positions  of  the  agents  will 
make  a  descending  chain  of  convex  sets. 

According  to  the  above  results,  those  agents  are  at  the  corner 
of  the  global  convex  hull,  RDD  algorithm  give  nonzero 
control  value.  The  convex  hidl  made  by  neighbors  of  agent 
at  time  instant  is  Co  ( )t  {z^kikip  where 

a  N..  It  has  been  proved  that  U  (.)  is  not  a  corner  point 
of  '^unless  z  z  ■■■□z  0. 

Theorem:  All  the  agents  will  converge  to  a  point  as  t 


Now,  ()  z\j$fy.\  ^(JnfoErges  to  P  P 

i  1 , . . . ,  m  1 .  Let  us  assume  by  slight  abuse  of  notation 
that  Vt()  is  r  in  the  sensing  region  of  X  .  We  know 

U(ti  r  iR(t  /)-  P(t ) 

From  the  above  equation  we  can  say  that  U();)  goes 
asymptotically  to  VX ,  which  is  a  corner  of  the  convex 
hull  because  V  is  a  corner  of  the  convex  hull.  According  to 
that  condition,  this  can  not  happen  unless  all  the  points  are  at 


point. 


Now  for  any  other  agents  two  conditions  can  arise 


Proof:  ObCtht;  convex  hull  generated  by  the 

points  04()0f  tb§\U<arfy  corner  Vtfi 

coincides  with  the  position  pt() of,  at  1$ ast,  one  agent 

a  at  time  instant  t .  Since  Co  t  (  )  Co  t  (  ) ,  then  it  is 

possible  to  define  the  convex  set  Co  as  the  closure  of  the 
intersection  of  all  the  Qfor  tCo  tt  0. 

Let,  V  be  a  corner  of  Co  .  There  exists  a  sequence  of  corners 
VtQk pf  Co  t  (),  which  converges  to  V  as 


Algorithm  Random  Random  LP  solved  Movements 

perturbation  flag  using  of  agents 

in  orientation  setting  agents 


A 

Y 

N 

All 

All 

flrsfikyf  is 

Y 

similar  to  that 

y  . 

given  in 

[9 

Active  . 
|.  So  trom  i 

Active 
t  we  can 

say  that  all  the  agents  will  converge  to  a  poilit.  The  second 
case  is  not  possible  as  the  initial  graph  is  connected. 

4  Simulation  Results  and  Discussions 


t  .  Again,  there  exists  a  value  j  { 1,  . . N}  ,  such 
that  p  Q,  |  (yt  for  infinite  number  of  values  of  /  . 

Let  us  assume,  without  loss  of  generality,  that  j  1.  Now 
let  us  select  only  those  values  t  of  t  with  /  0, 1 ,  ...  for 

which  pt[)u  Vt()  for  all  l  0. 

Since  the  number  of  neighbors  of  agent  a  is  non-  decreasing 
as  t  .  The  maximum  number  of  neighbors  should  be  N . 
As  the  number  of  neighbors  is  finite  then  we  can  discard  a 
finite  number  of  values  of  /  and  by  remembering  the  index 
variable  we  can  say  that  such  number  is  constant  (number  of 
neighbors  will  not  change).  As  the  set  {(),pt  I  0}  ,  where 
a  N  is  infinite  set  of  Jn 

points  and  is  a  compact  set  then  it  has  at  least  an 
accumulation  point  P  .  Then  it  is  possible  [9]  to  extract  a 

sequence  of  points  converging  to  P  . 


4.1  Simulations  Results 


In  the  first  set  of  simulations  we  have  three  agents.  Initially 
we  consider  ((1,1),45  ),  ((5, 4), 135  ),  and  ((2,6),  45)as  the 

initial  position  and  orientation  angle  of  the  three  agents. 
Using  the  perturbation  technique,  with  normal  distribution  for 
and  a  scaling  angle  of  120  ,  the  agents  converge  to  a 

point  after  a  few  iterations  (see  Figure  9).  The 
convergence  criterion  for  terminating  the  simulation  was 

when  the  value  of  r  became  less  that  210  . 


Let  us  select  only  those  values  of  (?|br^which  ptfy^ 

belongs  to  that  sequence  and  let  us  rename 

the  variable  /  by  r  .  This  implies  that  OcpMv^rges  to 

P  as  r 


Figure  6:  (a)  Trajectory  of  the  agents  (b)  Reduction  in 
r  with  iterations 


In  Figure  10,  four  snapshots  for  10  agents  have  been  shown 
in  two  dimensions.  Connectedness  of  the  agent  placement 
was  defined  as  follows;  two  agents  are  said  to  be  connected 
if  they  are  within  2  units  of  each  other.  Although  the  initial 


placement  of  the  agents  was  random,  we  ensured  that  the 
agents  were  connected,  which  implies  that  there  exists  at 
least  one  path  between  any  two  agents.  One  can  see  that  all 
the  agents  eventually  converge  to  a  point.  The  system 
converges  when  the  maximum  distance  between  the  agents 
along  the  global  reference  frame  is  less  than  the  decision 

domain  distance  ,  since  the  agents  would  converge  to  a 
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d 

single  point  in  the  very  next  step.  In  Figure  1 1  we  placed 
10  agents  randomly  in  three  dimensional  space.  One  can  see 
that  all  the  agents  converge  to  a  point. 
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(b)  (d) 

Figure  7:  10  agents  placed  randomly  converge  to  a  point  in  2D 
between  Ando  and  algorithm  in  two  dimensions 


4.2  Comparison  Results 

The  comparative  study  between  Ando’s  algorithm  and  the 
RDD  algorithm  was  carried  out  in  this  work.  In  Table  2,  the 
comparison  of  computational  time  to  convergence  is  given. 
The  results  show  that  RDD  algorithm  is  superior  in  terms  of 
computational  time.  The  average  was  taken  over  25  runs.  In 
Table  3comparison  of  the  number  of  iterations  to  converge  is 
given.  The  iteration  number  in  case  of  RDD  algorithm  is  a 
little  higher  than  Ando’s  algorithm.  The  reason  behind  this  is 
the  decision  domain  in  case  of  RDD  algorithm  is  smaller  than 
Ando’s  algorithm.  However,  this  is  easily  offset  by  the  fact 
that  RDD  takes  much  less  computational  time  per  iteration. 

In  Table  4,5  we  made  a  comparison 

between  Ando  and  RDD  algorithm  in  3D  space. 

Table  2:  Comparison  of  convergence  time 
between  Ando  and  RDD  algorithm  in  two 
dimension 


Time  to  Converge 


« 

Agents 

<“> 

- 

20 

0.5856 

0.0329 

17.8065 

50 

3.0667 

0.1758 

17.1544 

100 

,2.1486 

0.5439 

22.3360 

15 

40.7101 

1.7395 

23.4031 

0 

Table  3:  Comparison  of  number  of  iterations 
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Figure  8:  10  agents  placed  randomly  converge  to  a 
point 
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Table  4  :  Comparison  of  convergence  time  between 
Ando  and  RDD  algorithm  in  three  dimensions 


Time  to 

Converge 

Number 

Ando 

RDD 

Ando/RDD 

of  Agents 

(secs) 

(secs) 

Ratio 

20 

36.70 

0.03 

1079 

50 

543.72 

0.12 

4531 

100 

2319.75 

0.45 

5074 

150 

4655.62 

1.63 

2854 

ICIUS201C 
Nov  3-5,  2C 
Bali,  Indonf 


Table  5:  Comparison  of  number  of  iterations  between 
Ando  and  algorithm  in  three  dimensions 
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5  Conclusions 

We  presented  and  analyzed  a  rendezvous  algorithm 
considering  a  rectilinear  decision  domain.  The  computational 
complexity  of  RDD  algorithm  is  low  compared  to  the  well 
established  Ando’s  algorithm  in  the  literature.  The  RDD 
algorithm  is  simpler  in  terms  of  few  computations  needed  and 
in  relaxing  the  restriction  on  the  movement  of  the  agents  to  the 
centroid  of  the  rectangle. 


[10]  J.  Lin,  A.  S.  Morse,  and  B.  D.  O.  Anderson.  The 

nd 

multi-agent  rendezvous  problem.  In  Proceedings  of  the  42 
IEEE  Conference  on  Decision  and  Control,  pages  1508-1613, 
Maui,  USA,  2003 


Acknowledgement 

The  authors  gratefully  acknowledge  the  AOARD/AFOSR 
for  their  grant  to  this  project. 


References 

[1]  T.  Bretl,  Control  of  Many  Agents  Using  Few 
Instructions,  Third  Robotics  Science  and  System  Conference, 
Atlanta,  2007. 

[2]  K.  Das  and  D.  Ghose,  Positional  Consensus  in 
Multi-Agent  Systems  using  a  Broadcast  Control  Mechanism, 
American  Control  Conference,  St.  Louis,  Missouri,  USA, 
2009,  pp.  5731  -  5736. 

[3]  K.  Das  and  D.  Ghose,  Randomization  Mechanisms 
Based  Positional  Consensus  in  Homogeneous  Multi-Agent 
Systems,  II Sc  Centenary >  International  Conference  on 
Aerospace  Engineering,  Bangalore,  India,  2009,  pp  1327  - 
1335. 

[4]  K.  Das  and  D.  Ghose,  Multi-agent  Rendezvous 
Algorithm  with  Rectilinear  Decision  Domain,  FIRA,  2010, 
Bangalore. 

[5]  H.  Ando,  Y.  Oasa,  I.  Suzuki  and  M.  Yamashita, 
Distributed  Memoryless  Point  Convergence  for  Mobile 
Robots  with  Limited  Visibility,  IEEE  Transactions  on 
Robotics  and  Automation,  vol.  15,  pp  818  -  828,  1999. 


ISBN  978-S 


